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Abstract The explosion of a supernova releases almost instantaneously about 10®^ 
ergs of mechanic energy, changing irreversibly the physical and chemical proper¬ 
ties of large regions in the galaxies. The stellar ejecta, the nebula resulting from the 
powerful shock waves, and sometimes a compact stellar remnant, constitute a super¬ 
nova remnant (SNR). They can radiate their energy across the whole electromagnetic 
spectrum, but the great majority are radio sources. Almost 70 years after the first 
detection of radio emission coming from a SNR, great progress has been achieved 
in the comprehension of their physical characteristics and evolution. We review the 
present knowledge of different aspects of radio remnants, focusing on sources of the 
Milky Way and the Magellanic Clouds, where the SNRs can be spatially resolved. We 
present a brief overview of theoretical background, analyze morphology and polariza¬ 
tion properties, and review and critical discuss different methods applied to determine 
the radio spectrum and distances. The consequences of the interaction between the 
SNR shocks and the surrounding medium are examined, including the question of 
whether SNRs can trigger the formation of new stars. Cases of multispectral com¬ 
parison are presented. A section is devoted to reviewing recent results of radio SNRs 
in the Magellanic Clouds, with particular emphasis on the radio properties of SN 
1987A, an ideal laboratory to investigate dynamical evolution of an SNR in near real 
time. The review concludes with a summary of issues on radio SNRs that deserve fur- 
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ther study, and analyzing the prospects for future research with the latest generation 
radio telescopes. 

Keywords ISM; Supernova Remnants • radio continuum: ISM • radiation mecha¬ 
nisms: non-thermal • ISM: cosmic rays 


1 Introduction 


In 1919, the Swedish astronomer Knut Lundmark put forward the hypothesis that 
“along with ordinary novae, on rare occasions stars flare up that are tens of thousands 
of times as bright at maximum”. Zwicky and Baade proposed in 1934 that such stars 
be called supernovae, and though in the opinion of Shklovsky the name was rather 
absurd, it was rapidly popularized and now universally designates an event of stellar 
explosion ( Shklovsky |1978| l. 

Supernove (SNe) can be broadly classified into two big groups depending on 
the explosion mechanism; the Type la are the result of a runaway thermonuclear 
explosion of a degenerate carbon-oxygen stellar core (most likely a white dwarf in a 
binary system). The specific progenitor systems and the processes that lead to their 
ignition have not yet been clearly identified for these SNe as recently reviewed by 


Maoz et al. (20141. The other big group of SNe involves Types Ib, Ic and II, which 
are the product of gravitational core-collapse of massive stars (initial mass M> SM©) 
that have exhausted all their nuclear fuel. If the stellar core contains between 1.4 and 
3 Mq, the compact remnant is a neutron star, while if the mass of the collapsed core 
is larger than 3 Mq, a black hole is formed. As summarized by Smartt ( 2009| l there 
is a diversity of evolutionary scenarios, where metallicity, binarity, and rotation may 
play important roles in determining the end states, what kind of compact remnant 
leave, etc. As a general introduction to the radio SNRs and their connection with the 
explosion mechanisms that give birth to them. Fig. [T] summarizes the types of SNe, 
their precursors and possible remnants. 

Independently of the type of explosion, about 10®^ ergs of mechanic energy are 
suddenly deposited in the interstellar medium and several tens of solar masses of 
stellar material are ejected. The outer layers of the star blow off in all directions and 
the enormous explosion power imparts great velocities to the portions of the expelled 
envelope. A nebula is formed that expands at a speed that can reach 5000 to 10000 km 
s“^. This high expansion velocity is the main sign that distinguishes the remnants of 
supernova outbursts from other nebulosity. Hundreds or thousands of years later, the 
ejected material will begin to be slowed down by the ambient medium, the velocity 
will start to fall to hundreds or even tens of km s“^, and ultimately will disperse 
and merge with the surrounding gas. But through thousands of years a distinctive 
nebula, the supernova remnant (SNR) persists and can, eventually, be detected across 
the whole electromagnetic spectrum from radio to gamma-rays. The neutron stars, 
the compact objects left by most of the supernova explosions, will continue radiating 
energy for many more millions of years. 

Before the advent of radio astronomy only two SNRs were known, the Crab Neb¬ 
ula and Kepler’s SNR. Radio observations played a crucial role in the discovery and 
investigation of SNRs and their environs, revolutionizing the knowledge in the field. 
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Fig. 1 Diagram summarizing the types of SN explosion, their precursors and possible radio remnants. NS 
is the acronym of neutron star, BH of black hole and PWN of pulsar wind nebula. Mass limits are taken 
from|SmarttH2009y Mass quoted for NS and BH correspond to the compact core 


In 1948, at the very opening era of radio astronomy, the British astronomers Sir 
Martin Ryle and Graham Smith detected an unusually bright radio source, Cassiopeia 
A (Cas A). The flux of this source was comparable to the radio flux of the quiet Sun 
at meter wavelengths, and it seemed to come from a place where no optical emission 
was apparent. Faint, filamentary optical emission was only later discovered by Baade 
and Minkowski ( Baade & Minkowski||1954 l from plates acquired in 1951 with the 
200-inch Palomar telescope. But soon after 1948 several radio sources associated 
with SNRs were found in our Galaxy; the Crab Nebula, the remnants of the Tycho 
and Kepler supernovae, and a filamentary nebulae in the constellation of Cygnus. 
Then arose the fundamental question of the nature of those radio waves, since their 
radiation had nothing in common with thermal black body radiation, but rather was 
more akin to the nonthermal Galactic radiation detected in the early years of the 
radio astronomy by Karl Jansky and Grote Reber Pansky|1933[|Reber|1944| l. There 
was clearly a complete disagreement between the observed radio spectra of SNRs 
and those of thermal radio sources. The correct idea explaining radio emission of 
SNRs was proposed in 1950 by Alfven and Herlofson ( Alfven & Herlofson|| 1950) 1, 
and independently by Kiepenheuer ( |Kiepenheuer 19501. The answer was synchrotron 
radiation. The synchrotron nature was later confirmed by the detection of polarization 
in the Crab Nebula by Mayer et al. (1957 1 . 

In the following years, radio continuum surveys carried out at different frequen¬ 
cies were the main tool for identifying new SNRs through their non-thermal spec¬ 
trum, the “fingerprint” in this energy band. Whenever the instruments improve in 
sensitivity and angular resolution, new sources are found and this held is in continu¬ 
ous progress. From the very first lists of SNR candidates published by|Aizu & Tabara| 
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1967 1 and [Poveda & Woltjer| ( |T 968 | ) (with 25 SNRs listed), superseded by [Milne] 


1970 1 (who increased the list to 94 members in the Milky Way and 3 in the Large 


Magellanic Cloud), Downes ( 1971| l, [riovaisky & Lequeiix ( 1972| l, [Clark & Caswell 
( 1976| l, numerous works reported the discovery of new sources and refined their clas¬ 
sifications. A step forward came from the contribution of Brogan et al. ( 2006| l, who 
based on high-quality low-frequency radio observations disclosed the presence of 31 
new SNRs in the inner Galaxy, increasing in about 15% the total number of Galactic 
SNRs known by then. At the present time there are 2 94 firmly clas sified SNRs in 
our Galaxy, as compiled and permanently updated by Green ( 2014| P of which the 
large majority (~ 95 %) are radio sources. Also, Ferrand & Safi-H^b] (2012l keep 
an updated census of high energy observations of Galactic SNR^ listing their phys¬ 
ical properties along with a summary of observations of these remnants with various 
X-ray and 7— ray observatories. 

It has to be noted that after years of searching, the total number of detected SNRs 
in the Milky Way is significantly smaller than expected. A number at least 3 times 
greater of SNRs is statistically predicted on the basis of OB stars count, pulsar birth 
rates, SN rates in other Local Group galaxies, and predicted lifetime of radio syn¬ 
chrotron emitting sources in the sky. Such deficit is generally attributed to selection 
effects ( [Brogan et al.|2006) l, when old, faint, large remnants, as well as young, small 
sources, remain below the threshold in sensitivity and/or spatial resolution of the large 
Galactic surveys performed up to now. 


2 A brief overview of the theory needed to understand radio emission from 
SNRs 


Many good text books present the physical background of synchrotron emission in the 
radio astronomical context (e.g. ^^icholczyk|1970 Moffet 1975 Verschuur & Keller 


mann| |1988[ |Harwit|1988 | |Rohlfsjl990| Rohlfs & Wilson |1996 i, and more recently 

a good synthesis is presented in the review of SNRs at high-energies by Reynolds 
( 2008j l. Here we summarize just the basic results needed for the interpretation of ob¬ 
servational data with a brief guide to infer intrinsic properties and physical conditions 
of a source from the observed parameters. Readers not interested in the background 
aspects can jump to Section]^ 


2.1 Synchrotron radiation: 


A single relativistic electron moving in an external magnetic field B, will emit con¬ 
tinuum spectrum radiation. The frequency near which the emission reaches its maxi¬ 
mum intensity is called the critical frequency for synchrotron emission and is defined 
as: 


Vc 


3e 

47rmc \mc^ J 


CiB±E^ 


( 1 ) 


* http://www.mrao.cam.ac.uk/surveys/snrs/ 

^ http://www.physics.umanitoba.ca/snr/SNRcat/ 
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where the constant Ci = 6.266 x 10^® in cgs units or Ci = 16.08 when v is in MHz, 
B in microgauss and E in GeV ( Moffet|1975| l. For example, if we observe radiation 
at ~ 1500 MHz and assume a typical B field of ~ 10 /iG, we can conclude that such 
radiation was produced by electrons with an energy of 3 GeV. It can be shown that 
the time required for an electron to lose half of its energy is fi /2 = Cz/{B^ Eq) 
where C 3 = 8.35 x 10® yrs when B is in /rG and E in GeV. Then the electron of 
our example would radiate half of its energy in about 30 million years. For the Crab 
Nebula, with B^ 500 pG the time to lose a substantial portion of its energy for an 
electron radiating in radio (500 MHz) is ^ 100000 yrs, for an optical photon (600 nm) 
about 100 yrs, and for X-rays (4 keV) ^ 2.4 yrs. Such short lifetimes (much shorter 
than the age of the Crab Nebula) implies that the pulsar is permanently feeding energy 
into relativistic electrons. 

Actually relativistic particles are not alone but in an ensemble. The synchrotron 
radiation that an observer detects from a particular volume element of a radio source 
comes from all the electrons with the same pitch angle. We can assume that the 
source consists of a volume V containing a tangled magnetic field with average 
strength B, in which there are electrons with an energy density with power law dis¬ 
tribution (as the empirical evidence of the cosmic rays close to the Earth shows) 
n{E)dE = noE~^dE between a range of energies Ei < i? < i 72 to avoid diver¬ 
gences in the extremes. In this case it can be shown that the emitted radiation has 


a spectrum oc v 


(t-1) 


with a = 


(7-1) 


the emission spectral index . It 


can be easily concluded that the power law of the SNRs radio spectrum simply re¬ 
flects the power law energy spectrum of the relativistic particles responsible for the 
radio waves that we observe. After carrying out the corresponding calculations, it can 
be shown that the flux density of a synchrotron source with a volume V follows the 
expression; 

'6.26 X IQi®' 


S{v) = 0.0I7a(a)VS(“+i) 


v{Hz) 


Jy, 


( 2 ) 


where a{a) takes values like 0.283 for a spectral index a = 0 ; 0.103 for a = 0.5; 
or 0.085 for a = 0.75 ( Rohlfs||1990|l. The convention adopted here is that the flux 
density oc v~‘^ following Rohlfs (19901. 

During the evolution of a remnant the magnetic held strength will decrease pro¬ 
ducing a secular decay of the flux density for young radio sources ( |Shklovsky 1 1 960a| l . 
The exact rate depends on the held configuration, but assuming that the magnetic flux 
remains constant, B(R) is proportional to R~^ for an adiabatic expanding nebula of 
radius R. Then, the flux density should decrease with R as Si, oc R~®'>', which can be 
expressed as a function of time as Sjy oc Therefore, it is expected that the flux 

density decreases with time as: 


SJS, = -^. (3) 

bt 

This result is a rough approximation for a young source evolving into a uniform 
medium. A more rigurous result can be obtained by numerically computing the radio 
emission along with the hydrodynamic evolution of the SNR. In some cases for very 
young remnants (less than ~ 100 yrs), the radio emission can follow the opposite 
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behavior and increase the flux with time. This is the case for example of the youngest 
SNR observed in our Galaxy, Gl.9+0.3, whose flux density increased by a factor of 
1.25 over 13 yr at 1.4 GHz ( Green et al.|2008| l, and the SN 1987A in the Magellanic 
Cloud where the flux density increases in an exponential way with an average annual 
rate (measured at the year 20 since the SN explosion) of the order of 15% at different 
radio frequencies (Zanardo et al.|2010)l. 


2.2 Particle acceleration; 


An important question is what accelerates the particles to relativistic speeds in a SNR. 
If the remnant contains a pulsar, the answer is clear, it is the pulsar which supplies 
freshly accelerated electrons over the full lifetime of the SNR. But what happens with 
SNRs that do not have a neutron star in their interior? If the amount of material picked 
up by the supernova shock is large compared to the ejecta mass, then the particles can 
be accelerated in the shock waves. 

The basis of all theories to explain acceleration mechanisms in shock waves is that 
particles can gain energy in collisions with irregularities of the magnetic held. |Bell| 
( 1978a|b| l andjBlandford & Ostriker (19781 proposed that the most efficient process is 
diffusive shock acceleration (DSA), where electrons gain energy after multiple cross¬ 
ing through a shock wave ( Fermi|1949| l. For energetically unimportant test particles 
which do not influence the flow structure (the so-called test-particle limit), the syn¬ 
chrotron spectral index is fixed by the shock compression ratio r by a = 3/(2(r — 1)). 
In the case of strong shocks with a compression ratio of 4, this mechanism predicts 
a spectral index a = 0.5 for the radio flux density S. Coincidences and departures 
from this theoretical prediction are discussed in Section]^ For extensive discussions 
and reviews of standard and non-linear particle acceleration theories, see for exam¬ 
ple Blandford & Eichler| ( 1987 ; Reynolds (|2008| l; Maikov & O’C Drury P001| l; 
Reynolds (201 l|l; Jones (2011|l; Urosevic (2014|l, and references therein. 


2.3 Magnetic field; 


It was early theoretically suggested by |van der Laan| ( [l962 [) and later reinforced by 
Whiteoak & Gardner ( |I968 1 through a polarization study, that the magnetic field 
responsible for the synchrotron emission observed in SNRs was the interstellar field 
compressed by the explosion. The hypothesis was plausible because even a relatively 
weak compression of the few micro Gauss of the ambient magnetic field, would result 
in an observable radio source due to the strong dependence of volume emissivity on 
the magnetic field strength. 

Years later, however, indirect observational arguments, like the thickness of the 


X-ray rims in the SNRs Cas A, Kepler, Tycho, RCW 86 and SN 1006 (e.g. Jun & 
Norman|1996 Parizot et al.|2006| Ballet|2006 and references therein), and the rapid 
variability (on timescale of few years) observed in X-ray spots in G347.3-0.5 (SNR 
RXJ1713.7-3946) ( |Uchiyama et~ar]|2U07l l and in Cas A ( |Patnaude & Fesen|[2007l l, 
pointed to the fact that the magnetic fields must be from tens to hundreds of times 
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more intense than expected from adiabatic compression of the interstellar magnetic 
field, and some extra mechanism of amplification is required. Such mechanism must 
have little connection with the past of the SNR, since the seven studied cases have 
different progenitors, different explosion types (SNe types II and la), and contain or 
not a central neutron star, etc., but all of them require considerable magnetic field 
amplification to explain the observations. 

Several mechanisms have been proposed that amplify the magnetic field. |Schure| 
et al. (2012 1 and Reynolds et al. (|2012[) present reviews on observations and ' 


theories on magnetic field amplification in the presence of a cosmic ray population, 
showing that cosmic ray streaming can induce instabilities that act to amplify the 
magnetic field. Some of the proposed models include amplification due to flow in- 
stablities (of the Rayleigh-Taylor or Rankine-Hugoniot class) between the mean flow 
and clouds in the circumstellar and/or interstellar medium (e.g. Jun et al.|1995 Jun 


& Norm^|1996|l; amplification as a result of very efficient acceleration of nuclear 


cosmic 


rays at the outer shock (Volk et ar]|2005 and references therein); turbulent 


amplification driven by cosmic-ray pressure (Beresnyak et al. 2009 [ Drury & Downes 


|2012[ ); nonlinear DSA ( |Schure et al.|2012'l and references therein), etc. However, the 
global problem about the origin, properties and evolution of the magnetic field in 
SNRs, is still a matter of debate. 


2.4 The energy stored in particles and in magnetic fields: 


For an SNR, it is of interest to determine the energy content. This can be done by 
integrating over the electron spectrum between the energies Ei and i? 2 . If evidence 
of a cutoff is observed at one end or the other of the radio spectrum, the total electron 
energy can be expressed in terms of the cutoff frequencies i/i and 1/2 (the critical fre¬ 
quencies for the cutoff energies) and a, the spectral index of the source. By assuming 
that each electron radiates only at its critical frequency, we obtain an expression for 
the electron energy. 


Ue 


LCl^^ {2-2a) 
(1 - 2a) 


(if a ^ 0.5 or 1) 


(4) 


where L is the total luminosity of the source derived from the observed flux density 
S (L = Att(PS assuming isotropic radiation and that the distance d to the source is 
known), the constant Ci is as defined above, C 3 = 2.368 x 10“^ in cgs units and B is 
the average strength of the magnetic held. The magnetic held may not be measured 
directly, but its value may sometimes be inferred (at least an order of magnitude as 
is shown below). Observations at different frequencies provide an estimate of the 
spectral index a. 

The lower cutoff energy can be taken equal to the electron rest-mass energy (an 
electron that is not relativistic will not produce synchrotron radiation). Observations 
show departures from a power-law spectrum in the range from tens to few hundreds 
of MHz. Therefore a lower cutoff frequency of 10^ Hz is often assumed. For the 
upper cutoff, if the source under study has X-ray measurements, the break frequency 
derived from the intersection of radio and X-ray spectral slopes can be used as upper 
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frequency for radio emission. If there is no other indication, it is customary to use 
~ 10^° or 10^^ Hz as the upper cutoff. 

The energy of a source that emits synchrotron radiation can be mainly contained 
in two forms, as kinetic energy of the relativistic particles Upart and as energy stored 
in the magnetic field Umag (thermal radio radiation is negligible, therefore thermal 
energy will be small compared to the formers, and gravitational energy is probably 
the source of particles and magnetic energies, but almost nothing is known about it 
and is not considered). 

To estimate the energy stored in particles, we have to take into account that in 
addition to electrons there must be protons and other energetic baryons in the radiat¬ 
ing source (though heavy particles emit negligible amounts of radiation because they 
are accelerated much less by the Lorentz force). Therefore, Upart = where p 
is a factor that takes all other particles into account. From various methods, ry has 
been estimated to be 50 (the value from cosmic rays near Earth) or 100 (from models 
applied to radio sources in our Galaxy, e.g. |Burbidg^|1959| ). Fortunately, as shown 
below, this uncertainty does not have a strong effect on the fundamental conclusions. 

On the other hand, the energy stored in the magnetic field for a source of volume 
V is: Umag = VB"^ jSn. Therefore, 

Utot = Upart + Umag = TjALB + VB'^jSTT, (5) 


where A takes into account the constants Ci, C 3 and the shape factor in frequency. 

It can be noticed that the energy stored in particles, proportional to B~^U^ vvill 
dominate for small field strengths, while the magnetic energy (c>c B^) dominates for 
large fields. In consequence, the total energy Utot has a minimum for a magnetic field 
intensity 

( 6 ) 


for which Upart/Umag =4/3. Thus, the minimum total energy to produce the ob¬ 
served radio emission corresponds quite closely to an equipartition between relativis¬ 
tic particles and magnetic fields. If the magnetic field is exactly that of equipartition, 
then Umin = In this way, we can obtain a firm lower limit of the 

energy requirements of a synchrotron source. Examples of the use of these relations 
can be found in |Frail et al. ( 1996a| l for the PWN in W44, |Dubner et alI1 ( |2000| ) for the 
SNR W28, or in |Castelletti etaLjdSOOTl l for W44. 


The relations presented above are useful to get first estimates of energy content 
and magnetic field strength in radio SNRs, but it should be made clear that there is no 
physical justification for the energy components of the source being close to equipar¬ 
tition. It has been conjectured that motions in the plasma may stretch and tangle the 
magnetic fields and turbulent motions might also accelerate particles to high energies, 
thus taking the plasma closer to equipartition, but these are only conjectures and, in 
fact, these radio sources might be far from equipartition. 
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2.5 Dynamical evolution; 


Since the separation of the SNR evolution in four distinct phases originally proposed 
by Woltjer \\912\ , the scheme was maintained up to the present and is generally 
defined as follows: 

- Free expansion phase This phase occurs when the shock wave created by the 
explosion moves outwards into the interstellar gas at highly supersonic speed 
and compressed interstellar gas accumulates behind the strong shock front. This 
material is separated from the ejected stellar material by the contact discontinuity 
(a surface between two different materials with similar pressure and velocity but 
different density). Behind the contact discontinuity, a reverse shock starts to form 
in the ejected stellar material. After some time, the accumulated mass of the ISM 
compressed between the forward shock and the contact discontinuity equals the 
ejected mass of stellar material, and it starts affecting the expansion of the SNR, 
marking the beginning of; 

- The adiabatic expansion or Sedov-Taylor phase Here, the expansion is driven by 
the thermal pressure of the hot gas. In this stage, after the passage of the reverse 
shock, the interior of the SNR is so hot that the energy losses by radiation are very 
small. Since all atoms are ionized, there is no recombination, and the cooling of 
the gas is only due to the expansion. This stage has an exact self-similar solution 
( Taylor|1950[ Sedov|195^ . Later, as the SNR expands and cools adiabatically, it 
reaches a critical temperature of about 10® K, the ionized atoms start capturing 
free electrons and they can lose their excitation energy by radiation. The radiative 
losses of energy become significant, setting the end of the adiabatic expansion 
of the SNR. The efficient radiative cooling decreases the thermal pressure in the 
post-shock region and the expansion slows down. The SNR is entering in the 
so-called: 

- Snow plough or radiative phase when more and more interstellar gas is accu¬ 
mulated until the swept-up mass is much larger than the ejected stellar material. 
Finally, the shell breaks up into individual pieces, probably due to a Rayleigh- 
Taylor instability (hot thin gas pushes cool dense gas) and the SNR goes into the 
final phase. 

- Dispersion, as the expansion velocity decreases to values typical of the interstellar 
gas and the SNR disperses into the ISM. 


This scheme is an oversimplification and, as it was pointed out by Jones et al. 
( 1998| l distinct phases may be brief, may not occur at all, or may occur simultane¬ 
ously in different regions of the same remnant. A detailed description of the hydro- 
dynamical evolution of SNRs can be found in Lequeux ( 2005] l. 

|Chevalier| ( |1974| l developed a set of numerical models to describe the spherically 
symmetric expansion of a SNR in a uniform medium, providing semi-analytical ex¬ 
pressions to estimate expanding radius, age, shock velocity, etc. for different initial 
and ambient conditions. Although this work was later complemented with many dif¬ 
ferent variations to consider inhomogeneous surroundings, interaction with circum- 
stellar matter, etc. (e.g. |Chevalier|1982b[|Truelove & McKee|1999| l, it is still a classic 
that can be used as a first approach to derive physical parameters from observed 
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properties in SNRs. For example, for an initial energy Eq = 3 x 10^° ergs, ambient 
density no = 1 cm“^ and Bq = 3 x 10“® Gauss, the radius and shock velocity of 
an SNR at late times of evolution (after about 5 x 10"^ yrs) can be estimated within 
the model approximations using the following relations; i?n = 21.9 pc and 

^ ._n RQ , 1 ' ' 

Vsh 


= 66.5f^°®®kms-i, 


Vink 


(2012 I presents a good overview of different analyt¬ 


ical models developed for SNRs expanding in various environments. 


3 Radio morphology 


As already mentioned, the vast majority of SNRs in our Galaxy were first recognized 
from radio observations. In fact, out of 294 known remnants in our Galaxy, only 20 
have not been either detected in the radio band, or are poorly defined by current ra¬ 
dio observations ( Green|2014| l. The morphology and brightness distribution in SNRs 
contain important information about the nature of the SNR and its possible hydrody- 
namical evolution. 



Fig. 2 Detailed radio image of the Crab Nebula, obtained at 3 GHz with the Karl Jansky Very Large Array 
(NRAO) in the A-configuration. The image has a HPBW of x 0^^8 and the rms noise is better than 
30 /i.Jy (Dubner et al. 2015 in preparation) 
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Supernova remnants have been traditionally classified based on their radio mor¬ 


phology into three broad categories (see Weiler & Sramek 1988 for an historical per¬ 
spective on SNR naming conventions). They are: 


a) shell-type SNRs, whose appearance is characterized by a limb brightened shell or 
ring formed initially by the ejecta from the SN, and at later times also by swept up 
surrounding material. The diameter of the shell corresponds to the expanding shock 
wave produced by the explosion. In this case the particles responsible for the observed 
synchrotron radiation are believed to be accelerated at the shock front. 


b) filled-center or plerions, in which the radio brightness is centrally concentrated 
and is often linearly polarized. In this case, the accelerated particles and magnetic 
fields responsible for the synchrotron emission are injected by a pulsar created in the 
supernova event. This pulsar transfers the bulk of its rotational energy in a wind of 
relativistic particles and Poynting flux that interacts with the surrounding medium 
creating a synchrotron emitting nebula named “Pulsar Wind Nebula” (PWN). PWNe 
can emit non-thermal radiation over the whole electromagnetic spectrum. The Crab 
Nebula is a prime example of this class of remnants. Figure]^ shows a high-quality, 
sub-arcsec resolution radio image of the Crab Nebula at 3 GHz, obtained with the 
Karl Jansky Very Large Array (JVLA, NRAO) as part of a project of simultaneous 
multi-wave lengths investigation of this source ( Krassilchtchikov et al.|2014 Dubner 
et al. 2015, in preparation). See also Bietenholz et al. ( 2015| l for an image of Crab 
at 5 GHz with high resolution. It has to be noted that there are several plerionic 
SNRs, for which no pulsar has yet been found. For these cases the detection of the 
nebula is a strong indication that the powering source must be an undetected pulsar. 
Gaensler & Slane ( 2006 | l presented a complete review of PWNe as seen across the 
whole electromagnetic spectrum. 


c) composite, for which the SNR appears to have both a shell and an internal non- 
thermal pulsar driven nebula. As noted by Green ( 2014| l the term “composite” has 
been also used by some authors to describe SNRs with radio shell and centrally 
brightened thermal X-ray emission. Such SNRs are also known as “mixed-morphology” 
(M-M) SNRs (Rho & Petre 1998| l. Since the physical nature of the X-ray central 
component of these remnants is completely different from the non-thermal PWNe, 
we prefer to include them in a fourth morphological class. 


d) mixed-morphology (M-M), SNRs with synchrotron radio shell and central thermal 
X-ray emission. They have attracted more attention in the last years since several of 
them appear related to molecular clouds and to 7— ray sources, mostly of hadronic 
nature (when relativistic particles collide with dense ambient gas producing neutral 
pions, which decay into 7— ray ). The catalog by Ferrand & Safi-Harb ( 2012| l lists 
references for such associations. Vink ( 2012| l summarized their properties. Since new 
members have been reported in the last couple of years, to keep an updated census 
of these remnants we list in Tablel^the M-M SNRs in the Milky Way, together with 
their associations with molecular emission and 7— rays detections. 
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Table 1 Mixed-Morphology SNRs 


SNR 

Name 

Molecular 

OH (1720 MHz) GeV TeV Ref 



Material 

Masers 


GO.O-i-O.O 

Sgr A East 


Y 



1 

G6.4-0.1 

W28 

Y 

Y 

Y 

Y 

2,3,4, 5 

G31.9-HO.O 

3C391 

Y 

Y 

Y 


3, 6,7 

G33.6-H0.1 

Kes 79 

Y 


Y 


8,9 

G34.7-0.4 

W44 

Y 

Y 

Y 


10, 11, 12 

G41.1-0.3 

3C397 

Y 




13 

G43.3-0.2 

W49B 

Y 


Y 

Y 

14, 15, 16 

G49.2-0.7 

W51C 

Y 

Y 

Y 

Y 

17, 18, 19, 20 

G53.6-2.2 

3C400.2 






G65.3-H5.7 

G82.2-H5.0 

W63 






G89.0-H4.7 

HB 21 

Y 


Y 


21,22 

G93.7-0.2 

CTB 104A 






G116.9-1-0.2 

CTB 1 






G132.7-I-1.3 

G156.2-I-5.7 

HB 3 

Y 




23 

G160.9-I-2.6 

HB 9 






G166.0-I-4.3 

VRO 42.05.01 

Y 


Y 


24, 25 

G189.1-I-3.0 

G272.2-3.2 

IC443 

Y 

Y 

Y 

Y 

11,26, 27,28 

G290.1-0.8 

MSH 11-61A 






G304.6-I-0.1 

Kes 17 

Y 


Y 


29, 30,31 

G327.4-I-0.4 

Kes27 






G337.8-0.1 

Kes 41 

Y 

Y 



32, 33 

G344.7-0.1 


Y 


Y 

Y 

34, 35, 36 

G352.7-0.1 






37 

G357.1-0.1 

Tornado 

Y 

Y 

Y 


38, 39 

G359.1-0.5 


Y 

Y 



40,41 


References: (l)|Yusef-Zadeh et al.| jl996t ; (2)|Re ach et al j2005); (3 )|Frail et al.H1996b^; (4 ) |Aharonian 


et al.|j2008); ( 5)lHanabata et al.H2 014); (6)[Reach & Rho| 1999^; (7) |Castro & Slane 


et al.| 

2009 

et al.f 

2010 


S ; (17) 
; ( 21 ) 


; (9)|Auchettl^rtaljj20^; (1 0)|Seta^rtaUjV99y ; (1 l)jciausseneta^ 


Koo et al. )|Reichardt et al.fp^ ; (23) |Routledge et al.H199ll; (24) |Huan ~ 

Thaddeus]ji986); (25) Araya j2013^; (26)|Su et'S^j2014); (27)pckermann et al.H2013^; (28JAcciOT^ 


j20T0) ;(8) 


Giacani 


19^(12)1 Abdo 

; (1 6)|Abdo et al.| 


koo & MoonH1997) ; (18)[Green et al.l 1997t ; (19) |Abdo et 31.^20097; (20 ) |Feinstein et al.| 


et al.M2QQ9) ; (29)|Combi et al |i2010> ; (30 )|Gelfand et al H2013) ; (3 1)|Wu et al.||2Qi 1); (32) |Combi et al. 


j2QQ8); (^ ) |Koralesky et al.H1998); ( 34) |Giacani et al.H2Ql 1); (3 5) |Abdo et al.ff2013t; (36 ) |Aharonian 


et al.H200^; (37)|Giacani et ^.112009 


et al. 


; (38) Lazendic et al. H2004h (39)|Castro et al.H2013|; (40)|Hewitt 
<2008); (41)|Uchida et al.|fl^. 


In the recent compilation of Galactic SNRs by [Green (20141, 79% of remnants 
are classified as shell-type (including the mixed-morphology), 12 % as composite, 
and 5% as plerions. The rest of the remnants do not fit into any of the aforementioned 
conventional types, defying the simple spherically symmetric expansion solution. 


With regard to shell-type SNRs, multi-frequency observations of SNRs carried 
out in the last years with increasing resolution and sensitivity have demonstrated that 
the actual morphologies of these objects are highly complex and less than 20 % of 
the shell-types have an almost complete circular ring appearance. In effect, in radio 
waves SNRs exhibit an ample variety of shapes, such as blow-out: in which, part 
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of the shell appears to have expanded more rapidly than the rest. An example of 
this type is the SNR VRO 42.05.01 (G166.0+4.3) (Fig.^upper-right) that at first 


sight appears to be two different sources adjoining in the plane of the sky (Landecker 


et al.||1982]l, but that HI observations demonstrated that it is a single SNR breaking 


out from a warm medium into a tenuous interstellar cavity (Pineault et al. 1987| 
Landecker et al. 1989| l; barrel-shaped or bilateral: SNRs characterized by a clear 
axis of symmetry, low level of emission along this axis, and two bright limbs on 


either sides (Gaensler 1998 

1 . Typical examples of this type are G296.5+10.0 (IGiacani 

et al. 20001 and SN1006 ( 

Reynolds & Gilmore 1986| Petruk et al. 20091 (Fig. |3 


upper-left)', multi-shells: remnants with two or more overlapping rings of emission. 
Examples of this type are G357.7—0.1 (|Manchester||1987|l, 3C400.2 (|Dubner et al. 
[T994| ), G352.1-0.1 ( |Giacani et al.|2009| l (Fig.|3|hoftow-le/0. 

Another interesting morphology is created when a central compact object injects 
into a SNR a flow of particles along collimated bi-polar jets. This appears to be 
the origin of the distorted elongated shape observed in some SNRs, as in the SNR 
W50 (Dubner et al.|19^ Fig.^bottom-right), in Puppis A (Fig. 9 in Castelletti et al. 


20061, and in Cas A (Hwang et al. 20041. In the case of the W50/SS433 system. 


the radio observations confirmed the connection between the sub-arcsec relativistic 
jets from SS433 and the extended helical nebula over 5 orders of magnitude in scale 
( Dubner et al.||199^ . A similar origin was proposed for the “ears” observed in the 
SNR Puppis A, where the central neutron star would be responsible for the production 
of collimated outflows that impact on the shell. In the case of Cas A, the marked bi¬ 


polar asymmetry, in this case revealed in X-rays, was explained by Hwang et al. 


( 2004| l as the result of opposite symmetrical jets produced deep within the progenitor 
of the SN explosion. The nature of this morphology is the subject of observational 
and theoretical efforts to explain it. 

In the investigation of SNR morphologies, an additional complication is that the 
observed shape is a two-dimensional projection of a three-dimensional object and 
depends on their orientation with respect to the line of sight and projection effects. 
As a curiosity useful to illustrate this point, we reference the animation of a 3-D vision 
of the optical emission of the Veil Nebula SNR presented by J. P. Metsavainic]^ An 
example of the complications inherent to projection effects when studying the origin 
of the radio morphology is the case of the Galactic SNR G352.7-0.1 (Fig.^bottom 
left) whose radio emission projected in the sky plane looks like two concentric rings. 
This shape has been alternately proposed to originate in a “barrel-shaped” type of 
SNR ( |Giacani et al.|2009| l or from a blow-out scenario, where the SN explosion took 
place near the border of a molecular cloud (|Toledo-Roy et al.|2014|l. 


3.1 Can the SN type be constrained on the basis of the radio morphology of the 
SNR? 

This has been a challenging question for decades. The great diversity of shapes ob¬ 
served in radio SNRs, reflects not only different properties of the progenitor star and 
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http://astroanarchy.blogspot.com 
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Fig. 3 Radio images of SNRs with a variety of morphologies. Upper left: SN1006 at 1.4 GHz as taken 
from |Petmk et al.1j2009) ; Upper right: VRO 42.05.01 at 1.4 GHz as extracted from the CGPS public 
database; Bottom left: G352.7-0.1 at 5 GHz from|Giacani et al.H2009J; Bottom right: W50 at 1.4 GHz 
taken from |Dubner et al.H1998) 


of explosion mechanisms, but also echoes the properties of the ambient magnetic 
field and the matter distribution in the circumstellar and interstellar medium. Disen¬ 
tangling the different causes is a complex task that remains as a central question. 

Several attempts have been made to infer the type of supernovae on the basis 
of the observed remnants. In the X-ray domain, [Lopez et al. p009 20111 devel¬ 
oped an observational method to characterize the type of explosion of young SNRs 
by measuring global and local morphological properties of the X-ray line and ther¬ 
mal emission in numerous young SNRs in our Galaxy and in the Large Magellanic 
Cloud, finding that the remnants of Type la SNe have statistically more spherical and 
mirror-symmetric thermal X-ray emission than SNRs coming from core-collapse ori¬ 


gin. These studies were later extended to infrared (IR) morphology by Peters et al. 


( 2013| l through the investigation of the symmetry of the warm dust emission, conclud¬ 
ing again that Type la SNRs are statistically more circular and mirror symmetric than 
core-collapse SNRs. This is explained as due to different circumstellar environments 
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and explosion geometries of the progenitors. Additional ways of identifying the type 
of SN have been proposed using the X-ray emission, distinguishing la remnants from 
core-collapse ones by virtue of their ejecta composition ( [Hughes et al.|1995 1 (Fe-rich 
and O-poor SNRs are likely la, while SNRs dominated by O and Ne lines with weak 
Fe L emission are likely core-collapse ones). More recently, Yamaguchi et al. (20141 
presented an observational diagnostic to discriminate between progenitor types based 
on the Fe K-shell X-ray emission. The authors find that in remnants from SN la, the 
Fe-rich ejecta is significantly less ionized than in remnants from core-collapse SNe. 
Their results also indicate that there is a strong connection between the explosion 
type and the ambient medium density. 

The morphological criteria that seems to work for X-rays and IR to infer the 
SN type from the SNRs characteristics are clearly useless for radio SNRs. If we 
blindly apply the “circularity” criterion to Cas A and Tycho’s SNR, two circularly 
symmetric young radio remnants, both of them would be classified as SN la, but Cas 
A is core-collapse and Tycho’s is SN la. On the other hand, the SNRs SN 1006 and 
G296.5H-10.0 are two perfect examples of “mirror-symmetric” sources, but SN 1006 is 
a SN la, while G296.5H-10.0 has a central compact object suggesting a core-collapse 
origin ( Gaensler|1998| Harvey-Smith et al.|2010| l. 

One important reason that complicates the connection of a radio SNR with its 
precursor is that while X-rays may retain information about the characteristics of 
the exploded star, the complexity of the interaction between the shock front and the 
ejecta, circumstellar and interstellar matter, can soon mask this information in the 
radio emission. Once the shock front sweeps up a certain amount of ambient gas, 
the radio synchrotron emission ignores the explosion properties and it is mostly con¬ 
ditioned by inhomogeneities in the surrounding medium, hydrodynamic instabilities 
in the flow, turbulence behind the shock, effects of magnetic fields, etc. (e.g. e 


lier|1982a]b Dwarkadas|2005[). Sometimes, the traces left by mass loss episodes of 


the stellar progenitor can help to identify the class of supernova. The presence of a 
neutron star inside the remnant and/or the existence of the pulsar wind nebulae is an 
unquestionable evidence of a core-collapse event. The absence of it, however, does 
not prove anything because pulsars have high kick velocity (e.g. [Hansen & Phinney| 


19971 and can be outside the SNR far from the explosion site, also central compact 


objects can be radio silent, etc. An indirect indicator of core-collapse SN is hnding 
the SNR very close to or inmersed in a molecular cloud that might be the birthplace 
of a massive star that ended its life as SN Ib,c or II. 

In summary, the radio morphology alone does not provide a useful tool to distin¬ 
guish between different types of SNe. 


4 Polarization 

Radio polarization observations of SNRs provide essential information on the de¬ 
gree of order and orientation of the ordered component of the magnetic held, which 
themselves inhuence the morphology of the remnants and the intensity of radiation. 

As the radio emission in SNRs is primarily synchrotron, the radiation is linearly 
polarized; therefore from the observed polarization electric vector, the direction of 
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the orthogonally aligned magnetic field can in principle be determined. However, in 
practice, the observed polarization can be highly compromised by Faraday rotation 
since the electric field vector of the radiation is rotated during the propagation inside 
the SNR and the interstellar medium. In the simplest case, the angle of rotation ip 
varies proportionally to the square of the wavelength A, such that ip = RM A^, with 
the constant of proportionality, or the rotation measure RM, defined as 


RM(rad m ^) = 0.81 J N{cm (/rG')dl(pc) (7) 

where N is the thermal electron density and i3|| the magnetic field component 
along the line of sight 1, and the integral extends along the entire line of sight. The 
sign of RM is determined by whether Sy points toward or away from the observer. In 
general, observations at three or more wavelengths are necessary to measure the RM 
without ambiguity and determine the true position angles of the magnetic field lines. 

From the synchrotron theory, the intrinsic degree of linear polarization of the 
radiation emitted from electrons in a uniform magnetic field is independent of the 
frequency and given by P= (2a + 2)/(2a +10/3), where a is the radio spectral index 
of the radiation ( Ginzburg & Syrovatskii|19M||Pacholczyk|1970| l. For a non-thermal 
radio source with a spectral index a = 0.5, the fractional polarization can reach a 
maximum theoretical value of about 70%, but in practice a much lower polarization 
percentage is usually observed. 

The reduction of the observed degree of polarization arises from different physical 
and instrumental effects, namely beam depolarization as a result of variations of the 
rotation measure on spatial scales smaller than the antenna beam; bandwidth depo¬ 
larization, when the polarization angle changes across the receiver bandpass and the 
resulting non-parallel vectors are averaged; and depths depolarization (also known as 
differential Faraday rotation) when there is superposition of emission from different 
depths along the line of sight, either internal or external to the radiating source, which 
suffers different Faraday rotation. Under these conditions, the dependence of ip with 
A^ is no longer valid. To mitigate this effect several models have been proposed by, 
e.g. Burn (19661 and Sazonov ( 1973| l. For a discussion of basic depolarization mech¬ 
anism and instrumental effects, see also|Gardner & Whiteoak|(|1966|l,|Milne & Dickel| 
( fT975] l; |MiIne| ( fT980l l, and |R5di| ( |2006l l. 

Recently, a new tool, called RM Synthesis or spectropolarimetry, has been imple¬ 
mented to recover the polarization structure at multiple Faraday depths along a par¬ 
ticular line of sight with the additional advantage that minimizes n7r ambiguities and 


bandwidth depolarization (for a more detailed explanation see Brentjens & de Bruyn 


2005t[Heald et al.||2009| ). Its application has become viable thanks to technical and 
computational advances. Most application of this technique have been for interstellar 
or extragalactic studies (a good review can be found in Beck|201^ and its utilization 


for studying magnetic fields in SNRs is a very new field of research. Harvey-Smith 


et al. ( 2010| l applied this tool to the SNR G296.5H-10.0 using data taken with the 
Australia Telescope Compact Array at frequencies near 1.4 GHz. The Faraday rota¬ 
tion maps (Figure 1 in |Harvey-Smith et ar]|2010| l shows a highly ordered rotation 
measure structure with an anti-symmetric RM morphology across the remnant. The 
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authors proposed that the observed RM pattern is the imprint of an azimuthal mag¬ 
netic field in the stellar wind of the progenitor star. The expansion of the remnant into 
such a wind can account for the bilateral morphology of G296.5H-10.0 as observed in 
the radio and X-ray bands. 

Even if the radio observations of SNRs are carried out in optimum conditions to 
minimize the effects mentioned above, the observed degree of polarization in SNRs 
is still considerably lower than the maximum possible theoretical value. This is an 
indication that the magnetic fields are primarily disordered. In general, the polariza¬ 


tion degree has been found to be between 10 and 15 % (see references in Reynolds 


& Gilmorq|1993 i, with higher values between 35 and 60 % in some few exceptional 


cases, as : 


for example in some regions in the Vela SNR (|Mili^l9801l, DA 530 (|Lan-| 
decker et al.|19Wl l, G107.5-1.5 ( |Kothes|2003| ), and SN1006 ( [Reynoso et al.|2013^ 


At large spatial scale, it has been proposed that the intrinsic orientation of the 
magnetic field in SNRs, as inferred from the radio observations, shows a typical 
pattern depending on their age. The earliest observations of the young SNR Cas A 
showed a convincingly radial magnetic field with respect to the shock front (|Mayer 


& Hollinger|19681, while in the case of the old Vela SNR, the radio polarization map 


showed a near tangential direction in the brighter emission (Milne| 1968| l. Later on, 
polarization measurements carried out by Milne ( |1987| l over 27 SNRs confirmed that 
in young remnants the alignment of the magnetic field is predominantly in the radial 
direction, whereas the dominant orientation of the field in older remnants is parallel 
to the shock front or tangled. Subsequent observations have supported this picture 
( Landecker et al.| 19^91 Fiirst & Reich|2004t Wood et al.|2008[ ). In Figure]^ we show 
the intrinsic magnetic field distribution in the SNR Cas A (left) and in CTBl right), 
illustrating two extreme cases of magnetic field distribution predominantly radial and 
tangential, respectively. 

The general consensus on the tangentially ordered fields observed in older rem¬ 
nants is that they originate by compression in radiative shocks with large shock com¬ 
pression ratios. Regarding the radial component, however, the origin is still controver¬ 
sial. It has been generally attributed to stretching of Rayleigh-Taylor (R-T) instability 
fingers that occurs at the contact discontinuity which separates the ejecta from the 
shocked, more tenuous, circumstellar medium (CSM) (e.g., Jun & Norman 1996| l. 
However, polarization observations of SN1006 ( Reynoso et al.|2013| l show i^ial ori¬ 
entation near the forward shock (Fig. . This radio polarization study was carried 
out using data obtained with the VLA and ATCA instruments at 1.4 GHz with an 
angular resolution of about 10". The high quality of these data highlights the com¬ 
plex structure of the magnetic field distribution. A notable characteristic in SN1006 
is that even when the orientation of the magnetic field vectors across the SNR shell 
appears to be radial, a large fraction of the magnetic vectors lie parallel to the Galactic 
plane (whose orientation is perpendicular to the symmetry axis of the bright lobes). 
Reynoso et al. (|2013|l conclude from this evidence that the ambient magnetic field 


must be roughly parallel to the Galactic plane, and the SNR retains some knowledge 
of the orientation even after the passage of the shock front. Besides, while the de¬ 
gree of polarization in the two bright radio lobes of SN1006 is about 17%, a value as 
high as about 60% is found towards the faint SE border of the remnant. In brief, the 
brightest radio. X-ray and TeV features, the NE and SW lobes of SN1006, have the 
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Fig. 4 Radio polarization bars in B-field direction as obtained with the Effelsberg 100-m telescope. Left: 
SNR Cas A at 32 GHz; right: SNR CTBl at 10.55 GHz. (courtesy of W. Reich) 
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lowest polarization fractions (indicating the presence of a disordered, turbulent mag¬ 
netic field), while in the SE where the synchrotron emission is faint, the polarization 
is high (ordered field). In this case, the authors conclude that the most efficient par¬ 
ticle acceleration occurs for shocks in which the magnetic field direction and shock 
normal are quasi-parallel, while inefficient acceleration and little to no generation 
of magnetic turbulence are obtained for the quasi-perpendicular case, a result with 
implications for DSA theories. 

Over the past several years, a number of survey projects were launched in order 
to map the polarized emission from the Galaxy in great detail. Single-dish as well 
as synthesis telescopes located in the northern and southern hemispheres are being 
used. An updated list of the surveys of the Galactic polarized emission that have been 
conducted and being carried out, can be found in Tables 1 and 2, respectively, in 
|Landecker| ( |2012[ ), while a review of the radio polarization measurements from the 
beginnings to the present day was presented by Wielebinski ( 2012| l. As an important 
by-product, it can be mentioned that these surveys together with the total intensity are 
the main source for identifying new SNRs. Recently, 79 SNRs have been observed 
during the Sino-German polarization survey of the Galactic plane performed with 
the Urumqi 25-m telescope at 6 cm ( |Han et aT]|2014| l. Combining these data with 
observations made with the Effelsberg 100-m telescope at 21 and 11 cm, for the 
first time traced the magnetic field orientation for 23 out of 79 sources observed and 
identified two new SNRs, G178.2-4.2 and G25.1-2.3. 


5 Radio spectra 


The knowledge of the spectral behavior of SNRs is vital to understanding particle 
acceleration at the shock fronts and the role of SNRs as factories of Galactic cosmic 
rays. The study of the global and spatially resolved radio-continuum spectra of SNRs 
provides significant constraints on shock acceleration theories. As mentioned in Sec¬ 
tion 1^ in the case of strong shocks with a compression ratio of 4, DSA (diffusive 
shock acceleration) predicts for the radio flux density a spectral index a = 0.5. How¬ 
ever, when compared with the observed radio spectral indices, it is found that only 
50 out of the 294 SNRs listed in [Green (2014l’s catalogue have a between 0.5 and 
0.6, a number that only increases to 65 if the doubtful or poorly determined spectral 
indices are also considered. Eigurej^is an update of the summary of spectral indices 
presented by Reynolds ( 201 1[ ), after selecting the SNRs classified as shell-type and 
separating those sources with well known spectral index from those listed with a 
question mark, for which a simple power law is not adequate to describe their radio 
spectra. The members of the composite and mixed-morphology classes are included 
as shell-type radio SNRs. It is worth mentioning that among the SNRs classified as 
shell-type, almost 44% have poorly determined spectra, with a percentage as high as 
54% for sources located in the fourth Galactic quadrant, very likely due to observa¬ 
tional selection effects. These numbers, however, have to be considered with caution 
as the tabulated spectral indices in Green’s catalogue come from very diverse studies, 
with a wide range of (uncatalogued) errors. 
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Fig. 6 Histogram of spectral indices of shell-type Galactic SNRs as extracted from Green 


2014 . The solid 


bars correspond to firm spectral index determination, while the shadowed ones to SNRs with uncertain 
spectral indices 


Nevertheless, the fact that less than one fourth of the catalogued Galactic SNRs 
have the theoretically expected spectrum poses a problem to theorists. Different ex¬ 
planations have been proposed to address this. For instance, flatter spectra for increas¬ 
ing energy can be produced by shocks with low Mach number. The drawback of this 
hypothesis is that only a few SNRs would be expected to have such slow shocks, not 
enough to explain all the observed cases ( [Reynolds |2011| l. Also, second-order Fermi 
(stochastic) acceleration can play a role (Ostrowski 1999) and non-linear shock accel¬ 
eration (when accelerated particles influence the shock dynamics) can also produce 
flatter spectra ( Ellison & Reynolds |1 99 1[ Reynolds & Ellison|1992| l. 

Another important aspect to note is that theory predicts that the particle accel¬ 
eration must be very efficient in young SNRs, corresponding to free-expansion and 
especially early Sedov phases of evolution. This effect should translate in flatter spec¬ 
tra {a < 0.5) for young SNRs. However, observations show that young SNRs have 
steeper indices (e.g. a = 0.77 for Cas A, 0.58 for Tycho’s SNR, 0.6 for SN1006, and 
about 0.8, with extreme values in the range 1.1 to 0.3 for the remnant of SN1987A. 

In addition, the spectra of SNRs frequently show departures from a power law at 
the lower radio frequency extreme due to extrinsic or intrinsic processes that modify 
the electron energy distribution or the radio propagation. Three different physical pro¬ 
cesses that can affect the low radio frequency extreme have been proposed. They are: 
thermal absorption, synchrotron self-absorption, and the “Tsytovitch effect”. Ther¬ 
mal absorption occurs when the non-thermal SNR radiation traverses a region with 
thermal plasma, producing a low-frequency turnover in the spectrum. Synchrotron 
self-absorption takes place if the intensity of synchrotron radiation within the source 
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becomes sufficiently high, then reabsorption of the radiation through the synchrotron 
mechanism may become important, modifying the spectrum at low frequencies. The 
Tsytovitch effect can only be important in sources with very small size and weak 
magnetic fields (see Moffet [T975 for details). [Bell et alT ( 2011| l proposed that the 
spectrum can either steepen or flatten depending on the angle between the shock nor¬ 
mal and the large-scale upstream magnetic held. 

To narrow the gap between theory and observations and understand the physical 
meaning of the spectra, it is important to get good spectral index maps of many SNRs 
expanding in different environments and, if possible, in different evolutionary stages, 
to compare young and old population. The accurate radio spectral study of SNRs 
provides three different pieces of important information: (1) the global index, a key 
parameter to constraining particle acceleration theories as mentioned above, (2) the 
curvature of the spectrum, useful to test radiation mechanisms and also to separate 
intrinsic from extrinsic factors that may modify the spectrum, and (3) the existence 
of local variations within the remnant, that help to localize the sites where particle 
acceleration is taking place, the possible existence of radio PWNe (especially useful 
for the search of radio counterparts of PWNe discovered in other spectral regimes, 
like X- or 7- ray), and the presence of superimposed or embedded sources of thermal 
absorption. 

Since intrinsic or extrinsic spectral variations are often subtle, an important key 
is obtaining a large enough leverage arm in frequency space to tease them out. In 
the spectral studies it is essential to separate the contribution of three factors that 
usually overlap hiding the investigated properties: the intrinsic characteristics of the 
explosions, the contribution of the environment, and observational selection effects 
often imposed by the type of instrument used to acquire the radio data. This last one 
is not a trivial issue, since interferometric observations can over resolve the source, 
missing total flux information, but single dish observations that produce accurate to¬ 
tal flux measurements often confuse SNRs with nearby objects. An additional related 
issue is the problem of the proper subtraction of background emission, a serious issue 
since low-intensity emission extending over large angular distances (usually close to 
the Galactic plane) can result in a substantial contribution to the flux density of large 
SNRs. When analyzing the available radio spectra of SNRs, it is evident that most 
remnants are represented by only a few data points with substantial error bars, espe¬ 
cially at the lower frequencies. In what follows, we discuss the different techniques 
developed to overcome the various problems that can affect the spectral studies and 
some interesting results. 


Global spectrum To accurately determine global indices and their curvature, one 
technique is the T-T plot ( |Leahy & Roger|[l998| ), that fits pixel by pixel a linear 
relation to the brightness temperatures measured at two frequencies, deriving a tem¬ 
perature index from the slope (T^, = where /3 = a -I- 2). To investigate the 

possible curvatures, three or more frequencies must be used taking pairs of maps. 
This procedure guarantees that the derived spectral index, either total or local, is not 
sensitive to differences in the zero levels between the maps at the two considered 
frequencies. In any case, it is important that the images to be compared have similar 
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log (Frequency [MHz]) 


Fig. 7 Upper broadband spectral energy distribution from radio to 7 —ray showing the importance of 
the radio data to fit models in Tycho SNR {left, from|Giordano et al.|2012f and in SN1006 (right, from 
[Berezhko et al.|2Q12) . Bottom radio continuum spectrum for SNR 3C 391 showing a turnover in the 
spectrum at frequencies < 100 MHz, indicative of free-free absorption from thermal ionized gas along the 
line of sight (from[Brogan et al.|2005) 


angular resolution and sensitivity, whatever is the technique used to obtain the data 
(single-dish, interferometers or both combined). 

Kassim ( 1989|l carried out a pioneer work of observing numerous SNRs of the 


first quadrant at low radio frequencies and compiled existing observations, producing 
an atlas of radio SNR spectra, finding that about 10 out of the 32 observed SNRs show 
turnovers at low radio frequencies (below ^100 MHz), which was interpreted as due 
to the presence of a widespread, but inhomogeneous ionized absorbing medium along 
the line of sight, probably associated with the extended HII region envelopes (EHEs). 
Eigure [7] shows an example of spectrum with low frequency turnover (bottom), and 
two examples of broadband spectral energy distribution (SED) to illustrate the utility 
of the radio spectra to fit models to the 7— ray emission (upper). 


Spatial spectral variations The study in radio wavelengths of spatial spectral vari¬ 
ations inside a SNR can provide important information about the physical condi¬ 
tions of the interior plasma. One problem in the study of spectral variations across an 
SNR is the confusion with unrelated overlapping structures. To address this difficulty, 
|Katz-Stone & Rudnick|(|1997|l developed the tomographic technique, designed to iso- 
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late structures with different spectra. Spectral tomography involves making a gallery 
of maps using different test a, in which each tomographic image between frequen¬ 
cies and 1/2 is /((atest) = Ii — Features or regions which have a 

spectral index identical to the test value will vanish in the tomographic map. Spatial 
components that have different spectral indices will appear as positive or negative 
features depending upon whether the spectrum is steeper or flatter than the assumed 
test value. This method has been successfully applied to analyze spatial spectral vari¬ 
ations, for example in the SNRs of Tycho (Katz-Stone et al.|2000[), Kepler (DeLaney 


et al. 


et al. 


2002 1 , G292.0H-1.8 ( |Gaensler & Wallace||2003 1, and in Puppis A ( |Castelletti 


20061. The method is particularly useful to locate small-scale spectral varia¬ 


tions, providing a more accurate picture than the simple comparison between images 
at different frequencies. Figure|^shows the tomographic image obtained for Puppis A 
for atest =0.6, where it can be noticed a spectral pattern formed by short horizontal, 
almost parallel fringes with a alternatively steeper and flatter than the background, 
reproducing the “wave-like” morphology observed in the total power image along the 
NE, NW, and S borders of Puppis A. 


A significant breakthrough occurred with the first spatially resolved detection 
of thermal absorption towards the SNR 49B ( [Lacey et al.|200ll l, helping to explain 
the long mysterious presence of radio recombination lines towards this non-thermal 
source. 


Soon after, a number of additional cases of resolved thermal absorption came to 
light, as for example the observation of a cocoon of thermal material in the SNR 
3C391, marking the ionized boundary of the interaction of a SNR with a neighboring 
molecular cloud (Brogan et al.|2005]l. 


Another interesting example is the detailed spatial spectral study performed in 
Cas A, where Kassim et al. ( 1995| l showed that the data were consistent with absorp¬ 
tion by ionized gas inside the radio shell, probably related to unshocked ejecta still 
freely expanding within the boundaries of the reverse shock as delineated by X-ray 
observations. [DeLaney et alT] ( |2014| l extended these observations and their analysis to 
constrain the physical properties of a component of SNRs that has been heretofore 
very difficult to study. 


The spectral study carried out by [Castelletti et alT]p011a[ ) in IC443 revealed the 
existence of two different spectral components, both with flat spectrum (a < 0.25) 
but of distinct origin, one extrinsic and the other intrinsic. One of these components 
coincide with the brightest parts of the remnant along the eastern border and perfectly 
matches the region where ionic lines are detected in the J and H infrared bands. 
Such correspondence is the manifestation of the passage of a J-type shock across an 
interacting molecular cloud that dissociated the molecules and later ionized the gas. 
This ionized gas produces thermal absorption along the line of sight, resulting in the 
observed flattened spectrum. The other flat spectrum component is more fragmented 
and located near the center of the SNR, in spatial coincidence with a region with a 
molecular cloud and 7— ray emission (Fig. 0, suggesting that the origin in this case 
is the particle acceleration that takes place at the shock front. 
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Fig. 8 Tomographic image of Puppis A from 


Castelletti et al. 


'2006 


for Qtest 


a steeper than atest. while dark regions correspond to a flatter than atest 


= 0.6. Bright regions have 


6 Distances to SNRs 

A precise estimate of the distance to SNRs is essential to determine the physical 
parameters and understand their nature. It is, unfortunately, one of the most diffi¬ 
cult quantities to measure with accuracy. |Wilson| ( |1970| l expressed that “At present, 
distances to only a few SNRs are known with any reliability. Optical estimates are 
restricted to those SNRs which are rather near the Sun and, while nearly one hundred 
SNRs have been found (Milne 1970), distance estimates independent of any theory of 
SNR evolution are available for only a few sources”. Today, more than four decades 
later the situation have scarcely improved and only one-third of the catalogued Galac¬ 
tic SNRs has a distance estimate, several of them with large uncertainties or barely 
expressed as limits. 

The most precise determinations come from the combination of measurements of 
different independent quantities. |Trimbl^ ( |1973[ ) addressed the distance calculation 
for the Crab Nebula using 12 different lines of evidence through methods depending 
on the dynamics of the nebular expansion, on physical processes in the supernova or 
the remnant, and involving the properties of the interstellar medium between us and 
the nebula or the pulsar. This was an exceptionally well-studied case, and yet today 
it is affirmed that the true distance of Crab Nebula remains an open question owing 
to uncertainties in every method used to estimate its distance and the error amounts 
25% (d=2 ± 0.5 kpc). 

In general, for less famous sources, several methods have been used to mea¬ 
sure distances to Galactic SNRs, including the E—D (radio surface brightness-to- 
diameter) relation, kinematic distances through atomic and molecular absorption and 
emission, association with other objects with known distances (usually HII regions, 
OB associations, pulsars) assuming that they are neighbours. X-ray and optical obser¬ 
vations, etc. Each approach has its advantages and disadvantages, briefly discussed 
in what follows. 
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S—D relation Although the application of this method is very controversial, it is de¬ 
scribed for its historical value and because in some cases it may provide a hrst guess 
when no other distance indicator is available. It is only applicable to shell-type SNRs. 
The basic assumption is that since at a first approximation the mean radio sur¬ 
face brightness at frequency v, can be assumed as an intrinsic property of the SNR, 


it is a distance-independent observational parameter ( 

Shklovsky||1960a|b|l. The first 

theoretical X—D relation was derived by 

> 

o 

r/h 

1960a|) and the first empirical 

relations were derived by|Poveda & Woltjei 

' ( 1968|l. Clark & Caswell ( 1976|l describe 


that after reaching a maximum value shortly after the birth of the supernova, it may 
be expected that S decreases monotonically with time, while the outer linear diame¬ 
ter, D, of the expanding SNR will increase monotonically with time. The basic idea is 
then to construct an empirical relation using calibrators with distances known 
by other methods. It is expected that there is not too much scatter in this “evolution¬ 
ary track” . For the relation Ei, = AD^, the distance is derived as d cx eI,^^0~^ 
or, in terms of observable quantities, d cx sV^, where is the flux den¬ 
sity at the observing frequency v and Q is the angular diameter of the remnant. This 
formula also assumes that the SNR boundary is circular; when this is not the case 
it is commonly defined an equivalent diameter as 0 = 2-y/area/7r. As mentioned, 
departures from sphericity are not uncommon as they depend on environmental con¬ 
ditions and/or the orientation of the observer relative to the source. The ellipticity of 
typical remnants suggests that the variation of mean apparent angular diameter with 
orientation is sometimes as great as 10 %. 


The reliability of this method has been seriously objected, since there is no the¬ 
oretical basis for a single equation to cover the evolution of SNRs over their whole 
lifetime and expanding in a variety of ambient media (see discussions in Green|1984 


19911. Evolutionary paths may differ substantially from one SNR to another since 


they come from different stellar progenitors, probably experienced explosions with 
distinct physical characteristics, and may be expanding in interstellar media with a 
wide range of densities, from the evacuated interior of a wind-blown bubble, to a 
dense molecular cloud. However, in spite of the serious objections, this method has 
been widely applied for years, and there have been several attempts to improve it 
in different ways. The distances of the sources used as calibrators were refined by 
changing the kinematical estimates through the application of a more modern rotation 
curve for the Milky Way, by taking new data from associations of SNRs with pulsars 
or molecular clouds, or by comparing with new X-ray and optical data. The different 


estimated power law index /3 varies in the range -2.38 to -5.2 (see Case & Bhat- 


tacharya|1998 for a summary and comments on the different calibrations). Pavlovic 


et al. (2014) presented the most recent fitting based on a revised calibration sample 
consisting of 65 Galactic SNRs, deriving f3 = —5.2. Arbutina et al. ( 2004) showed 
that the only case were a single E—D correlation could be established was for the 
SNRs in the starburst galaxy M82, and suggested that a similar behavior could exist 
for a sample of Galactic SNRs associated with large molecular clouds. They con¬ 
clude that the density of the environment is probably a crucial parameter to regulate 
the diameter during the SNR evolution. 
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Kinematical method Basically, this method is based in the construction of an absorp¬ 
tion HI spectrum by subtracting an average spectrum obtained from an area projected 
against a region with strong continuum emission of the target SNR, from a spectrum, 
or average spectrum, of an adjacent background region. Using the radial velocity of 
the absorption peak and by applying a Galactic circular rotation modeQto convert ra¬ 
dial velocities into distances, it can be set a distance, or range of acceptable distances, 
for the SNR. Complications arise when the HI emission is patchy and cause spurious 
absorption features or when the continuum is faint and it is not possible to construct 


acceptable absorption spectra, with peaks noticeable at least at a 5-a level. Tian & 


Leahy (|2008 1 apply a combination of HI absorption plus CO emission to overcome 


some of the basic problems. Kinematical distances can also be derived from the study 
of the ISM around the SNR. When there is evidence of interaction between the SNR 
and the surrounding atomic or molecular gas (Section]^, the radial velocity at which 
the best signature is identified can be used to establish an approximate distance to the 
SNR. 


However, even when the radial velocity of the absorption or emission associated 
feature is clearly defined, there still exists the problem of the ambiguity of the Galac¬ 
tic rotation curve within the Solar circle, i.e. for each radial velocity there are two 
corresponding distances equally spaced on either side of the tangent point. Moreover, 
although with the help of other indicators a single distance can be established, the cir¬ 
cular rotation model still adds large intrinsic uncertainties (the model holds only for 
low Galactic latitudes and ignores systematic stream motions in the Galaxy, rolling 
motions in the Galactic arms, random cloud velocities, etc.) and the distances de¬ 
rived using kinematical methods have intrinsic (inherent to the method) uncertainties 
greater than ^ 25% - 30%. 


Distances derived from X-ray observations Kassim et al. (19941 proposed a formula 
for shell-type SNRs to estimate distances in the cases where it can be assumed that the 
SNR shell is in the adiabatic expansion phase, that the measured X-ray temperature 
gives a reliable estimate of the SN shock velocity, and that an initial energy of the 
order of Eq 10®^ erg is valid for all SNe. In these cases, the distance to a SNR 
could be derived as a function of the initial energy, the observed angular diameter 
of the SNR shell, the measured X-ray flux corrected for interstellar absorption, the 
thermal temperature of the X-ray emitting gas, plus a function that describes the 
power emitted by hot electrons in a low-density plasma via free-free emission (that 
depends of both the energy band of the emission and the temperature of the plasma) 
as Ds = 8.7 X 10®e[i'4p(Z\£;, The method was applied to a 

sample of SNRs detected with ROSAT, with the major uncertainty related to the basic 
physical assumptions rather than to measurement errors. 

Other application of X-ray data to constrain distances to SNRs was for example 
the one used in the study of the SNR RX J1713.7-3946 (G347.3-0.5) first discovered 
by its X-ray emission ( |Cassam-Chenai et al. 2004| l. The procedure consists of plotting 
the cumulative absorbing column calculated from atomic and molecular gas observa¬ 
tions as a function of the radial velocity and comparing with the absorbing column 


4 


Updated models for the Milky Way ai‘e presented in Bhattacharjee et al. (2014) 
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NH as derived from a fit to the X-ray data. The radial velocity at which these quan¬ 
tities are equal is later translated to distance through a kinematic model of Galactic 
rotation. 

In summary, the distances estimated for SNRs have to be considered with cau¬ 
tion. Even in the best cases when more than one independent distance determination 
can be used, the derived distances are still imprecise, either because of observational 
inaccuracies or because of doubtful assumptions involved in the formulae used to 
calculate them. 


7 Interaction of SNRs with the surrounding ISM 


The interaction between the strong SN blast wave and the surrounding interstellar 
medium has profound consequences on the remnant as well as on the gaseous inter¬ 
stellar matter. SN explosions are the main way of chemical enrichment of the ISM, 
while the distribution and physical conditions of the surrounding gas represent the 
primary physical constraint to the expansion of the SN shock. The investigation of 
the SNR/ISM interaction is not only necessary to improve our knowledge of SNRs, 
but also to understand the response of the interstellar gas to the enormous injection of 
energy and momentum that a SN explosion represents. Additionally, as already men¬ 
tioned in Sec.|^ the identification of physically associated interstellar gas serves to 
calculate the kinematical distance to the remnant. The study is also very important to 
understand the production of 7 —rays, and hence the role of SNRs powering Galactic 
cosmic rays, as clouds illuminated by the protons accelerated in a nearby SNR could 
be bright 7 —ray sources (e.g. Aharonian & Atoyan 1996 Berezhko & Volk||2000 


IGabici et al.|2009l|Butt|200^|Ellison & Bykov|201l|lT 

The ISM is formed by various components in different physical conditions. It 
is arranged in a variety of structures and phases, including big complexes of dense 
clouds, hot bubbles, sheets, walls, filaments, etc. (e.g. Lequeux 2005 [ Cox 2005) 1. 
Some of these structures are directly created by the SNe and SNRs. The environmen¬ 
tal characteristics play a crucial role in the shape, energetics, temporal evolution and 
destiny of SNRs. Before analyzing the consequences of the mutual interaction be¬ 
tween SNRs and the ISM, in what follows we briefly review the basic characteristics 
of the ISM in our Galaxy. 


Phases of the interstellar medium: 

- The molecular medium (MM), characterized by cold dense molecular clouds which 
are mostly gravitationally bound, with typical temperatures < 100 K, volume den¬ 
sities > 10^ cm“^, and volume filling factor f < 1%. A mass of about 1.5 x 10® 
Mq of molecular gas is distributed predominantly along the spiral arms of the 
Milky Way, occupying only a very small fraction of the ISM volume and within 
a narrow midplane with a scale height Z ^ 50 to 75 pc. 

- The cold neutral medium (CNM), distributed in rather dense filaments or sheets, 
with typical temperatures of ^ 100 K, volume densities ~ 20 - 25 cm“®, and vol¬ 
ume filling factor f ^ 2 to 4%. This phase is most readily traced by HI measured 
in absorption. 
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- The warm neutral medium (WNM), which provides the bulk of the HI seen in 
emission, with typical temperatures > 6000 K, volume densities ~ 0.3 cm“^ and 
volume filling factor f > 30%. 

- The warm ionized medium (WIM) with T ^ 8000 K, n ^ 0.03 - 0.3 cm“^, and f 
> 15%, ionized gas associated with HII regions, but also diffuse filling a consid¬ 
erable fraction of the ISM. 

- The hot ionized medium (HIM) with 10® K, n~ 10“® cm“® and f < 50%. This 
phase of hot gas is produced by supernova explosions, has a long cooling time and 
consequently a large fraction of the ISM is filled with this “coronal” gas. It can be 
viewed as the accumulated superposition of dissipated SNRs interiors integrated 
over the lifetime of the Galaxy (McKee & Ostriker|1977[). 


The cold, warm and hot phases are in global pressure equilibrium, while the 
molecular material is mostly confined to clouds which are held together by gravi¬ 
tation. The filling factor for each of the phases is highly uncertain, as is the topology 
of the ISM. In general, the ISM has structures on all scale lengths, from smaller than 
1 pc to larger than 1000 pc. Through HI observations of the Galaxy it has been known 
that the ISM is pervaded by holes surrounded by shells and supershells (|Heiles 1984 


IMcClure-Griffiths et al.|2002i|Suad et al.|2014D , as the result of individual or collec¬ 
tive action of stellar winds of massive stars and supernova explosions. The interior of 
these bubbles is filled with HIM. 

Molecular gas, on its side, is mostly (~ 90%) assembled in massive structures 
distributed in large clumps (GMCs, the giant molecular clumps with masses of ~ 
10"* — 10® Mq, diameters ^ 50 pc and average densities ^ 100 — 300 cm“®), 
clouds (M ~ 10^ Mq, diam ^ 5 pc, 300 cm“®), and condensed cores (possi¬ 

ble birth sites of new stars, with M ~ 10® Mq, diam ^ 2 pc, ^ 10® cm“®) (e.g. 


Williams et al.|2000| l. They are principally compose of by molecular hydrogen and 
dust and, in very little quantities, by some of the ~ 140 molecular species that have 
been identified in interstellar or circumstellar gas. 


7.1 Consequences of the interaction; 


As expected, the expanding SNR shock front undergoes a different evolution accord¬ 
ing to the phase of the ISM that it encounters in its expansion, and the dynamical 
evolution of a SNR is modified with respect to the simple description presented in 
Section Given a constant pressure in the SNR and that the postshock pressure 
scales with (where po is the ambient density and Vs the shock velocity) it is 
expected that the SN rapidly expands in the lowest density phase, while dense clouds 
are slowly crushed by lower velocity shocks. The energy can be conducted from the 
hot gas filling the SNR to embedded clouds, leading to their evaporation. A complete 
treatment of the expressions that govern the expansion of a SNR in these cases can 


be found in Tielens (20051. In the end, the exchange of mass, energy, and momentum 
between the different phases govern the dynamical evolution of the SNR, and the 
appearance of the SNR in the different spectral regimes reflects the structure of the 
surrounding medium. 
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Due to the short lifetimes of massive stars, most core-collapse SNe (resulting 
from explosions of types SNII, SNIb and SNIc) are located close to the molecu¬ 
lar concentrations where the precursors were born. Therefore a large percentage (as 
high as ^ 75%) of the Galactic SNRs are expected to be interacting with MCs. One 
important consequence of these interactions is that because of the compression and 
high temperatures of the SNR shocks propagating inside MCs, some chemical reac¬ 
tions otherwise not possible, can occur, creating new molecular species (e.g. |Tielens| 
20051. For example, a recent work by Dumas et al. ( 2014)1 reported the detection of 


SiO emission triggered by the passage of the W51C SNR shock. 

Slane et aL](2015 i describe the X-ray and 7 —ray signatures of the interaction of 


SNRs with molecular clouds and summarize the different criteria that can be used 
to establish the existence of a physical relation between an SNR and a cloud seen in 
projection against the remnant. To unambiguously establish whether an SNR is phys¬ 
ically associated with an interstellar cloud, removing confusion introduced by unre¬ 
lated gas along the line of sight, is not trivial and usually requires of several distinct 
criteria to demonstrate physical interaction. Basically, morphological traces along the 
periphery of the SNRs, such as arcs of gas surrounding parts of the SNR or inden¬ 
tations in the SNR outer border encircling dense gas concentrations. Usually, such 
features are indicating that a dense external cloud is disturbing an otherwise spheri¬ 
cally symmetric shock expansion. These initial signatures need to be confirmed with 
more convincing, though more rare, features like broadenings, wings or asymmetries 
in the molecular lines spectra ( Frail & Mitchell) 1998 1 , high ratios between molecular 
lines of different excitation state ( )Seta et al.| 1998) 1, detection of near infrared H 2 or 
[Fe II] lines (e.g. Reach et al.))2005)l, peculiar infrared colors (e.g. Castelletti et al. 


201 la)l, and the presence of OH (1720 MHz) masers. These masers were originally 


detected by Goss & Robinson ( 1968)1 in some SNRs, but their importance as proof 


of SNR/MC interaction was only recognized almost 30 years later (Claussen et al. 


)1997) [Frail & Mitchell|199^ , and their detection became the most powerful tool to 
diagnose SNR/MC interaction. In addition to be an ideal tracer of interaction, ob¬ 
servations of OH (1720 MHz) masers provide an accurate estimate of the magnetic 


field intensity of the post-shock gas in SNRs via Zeeman splitting (e.g. Brogan et al. 

[ 2000 ) 1 . 


Based on a combination of different techniques, Jiang catalogued a list of ~ 70 
Galactic SNRs candidate to be physically interacting with neighbouring MC^ of 
which 34 cases are confirmed on the basis of simultaneous fulfillment of various 
criteria, 11 are probable, and 25 are classified as possible and deserve more studies 
( Chen & Jiang|201^ . 

The basic method of investigating cases of SNR/ISM interaction is to survey the 
interstellar medium in a field around the SNR using different spectral lines, from the 
cold, atomic hydrogen emitting at A21 cm to the dense shielded regions of molec¬ 
ular hydrogen emitting in the millimetric and infrared ranges. The molecular gas is 
frequently studied through CO observations (usually the transition ^^C^®0 J: 1-0 at 
A 2.6 mm), since this molecule radiates much more efficiently than the abundant 
H 2 and can be detected easily. Later, the observed CO intensity has to be converted 


5 


http://astronomy.nju.edu.cn/~ygchen/others/bjiang/interSNR6.htm 
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into total H 2 gas mass. According to the standard methodology, a simple relation¬ 
ship can be established between the observed CO intensity and the column density of 
molecular gas, such that N(H 2 ) = XCO x W(CO J; 1-0), with the column density, 
N(H 2 ) in cm“^ and the integrated line intensity W(CO) in K km s“^. Bolatto et al. 


( 2013| l present a complete review of the theory, techniques, and results of efforts to 
estimate XCO, the “conversion factor” from CO into H 2 in different environments. 
In the Milky Way disk, the representative XCO values varies between 0.7 and 2.8 

(K km s“^)“^ depending on the technique applied for the estimate. A 
20 


X10^° cm ^ 


value XCO = 2 x 10^*^ cm"^ (K km s 




by Bolatto et al. ( 2013| l. 


) ^ with 30% uncertainty is recommended 


In HI (absorption and emission methods), numerous studies around SNRs have 
been carried out using single-dish and interferometric radiotelescopes, looking for 
traces of disturbances in the interstellar gas distribution caused by the SNR or its 


stellar progenitor (e.g. Dubner et al.|19981 Reynoso et al.|1999 Giacani et al.|2000| 

Velazquez et al.|2002 

[Koo & Kang|2004|[Paron et al.|2006 Leahy & Tian|2008 Lee 

et al.|2008 Park et al 

.|2013|l. An interesting approach using atomic gas studies was 


recently used in the historic remnant of SN 1006 by Miceli et al. (2014 1 , where an 


important connection between shock-cloud interaction and particle acceleration was 
demonstrated based on the comparison of X-ray with HI data. 


Many dedicated studies were conducted towards several Galactic SNRs in dif- 


ferent molecular transitions (e.g.lReach & Rho|1999[ 

Moriguchi et al.|20011 Dubner 

et al. 

2004 Moriguchi et al.|2005 Reach et al.|2005 

Jiang et al.|2010 Paron et al. 

2012 

Li&Chen|2012 Hayakawa et al.|2012 Fukui et al.|2012 Gelfand et al.|2013 

Kilpatrick et al.|2014 etc.). Of particular importance have been the molecular gas in- 


vestigations around the SNRs IC443, W28 and W44, which have been the subject 
of many observational studies. These were the first SNRs where the OH(1720 MHz) 
masers were observed. In the recent years, these sources became especially notorious 
as they are good examples of Galactic 7 —ray sources detected in the TeV range as 
seen by H.E.S.S. and other Cerenkov telescopes, and in the GeV range by the Fermi 
and AGILE satellites. 


The SNR IC443 (G189.1-1-3.0), because of its location in a relatively confusion 
free region of the outer Galaxy, is a text book case to analyze shock chemistry, and 
as such it has been thoroughly studied in many molecular transitions. From the first 
identification by [Cornett et al.| ( |1977| ), tens of works investigated the chemical and 
physical transformations introduced by the strong SNR shocks on the surrounding 
molecular clouds (see [Kilpatrick et al.||2014 for an updated summary of molecular 
studies towards IC443). The morphology consisting of two semi-circular shells with 
different radii is an indication of expansion in an environment with a marked density 
contrast. 1C443 is also an excellent example where the interaction of the SNR with 
a molecular cloud probably gave origin to the 7 — ray emission through a hadronic 
mechanism. In effect, an excellent concordance was demonstrated by[Castelletti et al. 


([201 lab by confronting the VERITAS very high energy 7 —ray radiation (Acciari 
et al.[[2009| 


with the J=l-0 cloud ( Zhang et al.[[20ld i (Fig. [^, and the region 
where the radio-continuum emission of IC443 shows flat spectral indices not caused 
by thermal absorption (see Section]^. 
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Fig. 9 The SNR IC443: left: radio-continuum at 330 MHz (in greys, from Castelletti et al. 2011a with 
7 —ray emission from VERITAS (contours), middle: VERITAS image of very-high-energy (VHE, E > 
100 GeV) 7 -ray emission jAcciari et al.|200^ ; right: Molecular CO emission (in greys, from |Zhang et al.| 
| 2010 |with 7 —ray emission overlapped in contours) 


The SNR W44 (G34.7-0.4) (Fig. [T0)l is another Galactic remnant whose interac¬ 
tion with a cloudy ambient medium has been deeply investigated, at first motivated 
by its morphology of a rather distorted shell with a significant flattening along the 
eastern edge, suggestive of an encounter with dense ambient medium. Since the first 
HI study by [Sato ( 1974|l, reporting the presence of a dense cold cloud in coinci¬ 


dence with W44, many observations have been performed in different molecular and 


IR lines (e.g. 

Dickel et al.|1976 

|Wootten||1977 

Seta et al.||2004 Paron et al.|2009 

jSashida et al. 

2013 Anderl et a 

.|2014|and references therein). The IR data showed 


that this SNR expands in a dense medium with n ~ 100 cm 


-3 


Giuliani et al.|( 201 1[) used constraints set by CO data from the NANTEN Observatory, 


the radio spectrum as obtained by |Castelletti et al.] ( |2007| l and optical data reported by 
IGiacani et al.| ( |1997| l, to demonstrate that the 7— ray emission detected with AGILE 
in the energy range 400 MeV — 3 GeV, was consistent with hadron-dominated mod¬ 
els. The possibility of star formation in the region of strong interaction between the 
expanding SNR shock and the adjacent molecular cloud (where there are IR sources 
with spectral characteristics of young stellar objects), is discussed below in Section 
7.2. 

Another interesting observational example of interaction can be found in the SNR 


Puppis A (G260.4-3.4), where CO and HI observations carried out by Dubner & 
Arnal ( 1988|l first revealed that the SNR is expanding in an inhomogeneous environ 


ment, with clouds towards the east and northeast borders. Later, [Reynoso et ah ( 1995| ) 
explored this SNR in HI with high angular resolution using the VLA (Eig. 11 Left) 


revealing the distribution of a pre-existing interstellar HI cloud closely following the 
borders of the SNR shell as seen in X-rays. Inside the SNR, Hwang et al. ( 2005| l in¬ 
vestigated through X-ray images and spectra the most prominent spots of cloud-shock 
interaction (mainly the bright eastern knot, known as the BEK), revealing what has 
been the first X-ray identified example of a cloud-shock interaction in an advanced 
phase. Erom new CO observations, Paron et al. ( [2008 1 showed evidence of what 
has been left of the molecular clump engulfed by the SN shock and that currently 
evaporates emitting in X-rays. Another interesting conclusion is that to the eastern 
boundary where the atomic gas is in touch with the SNR, the molecular emission is 
detached, indicating that the precursor radiation has dissociated the molecules of the 
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Fig. 10 Left: Radio continuum emission of the SNR W44 at 324 MHz (VLA image, in blue) with the 
infrared emission in the SNR region from Spitzer data at 8 pm (in green) and at 24 pm (in red) (from 
|Castelletti et al.|2Q07^ . The yellow spot in the southeast con'esponds to the HII region G034.8-0.7, an area 
where there is active star formation. Right: J = 1-0 distribution showing the giant molecular cloud 

GMC G34.8-0.6 interacting with the SNR W44. The cross shows the position of the infrared source IRAS 
18542+0114. The blue contours trace the silhouette of SNR W44 displayed in Galactic coordinates (from 
|Paron et al.|2009) 


adjacent cloud. It is also noticeable from Figure [TT]f?/g/if that in X-rays Puppis A has 
a “cellular” hlamentary structure with a honeycomb appearance, conhrming that it 
expands in a rich and complex environment. In addition, Figure showing 

the location of the radio border running outside of the sharp X-ray and IR limbs, is 
clear evidence that the radio radiation is the best indicator of the true position of the 
shock front. 



Fig. 11 Left: Distribution of the HI emission (in contours) around the SNR Puppis A, overlapping the 
ROSAT X-ray emission (in greys) as taken fromjPetre et al.|jl982|. Image extracted from|Reynoso et al.| 
Right: Two color image of Puppis A comparing the Chandra and XMM-Newton X-ray emission 
in the 0.7-1.0 keV band (in green) jPubner et al.|20f^ , with IR emission as observed with Spitzer at 
24pm (in red) ^Arendt et al.|2Q10) . The white contours show radio-continuum emission at 1.4 GHz (from 
ICastelletti et al'|2006t ' ' 
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7.2 Can SNRs trigger the formation of new stars? 


Since |Opik| ( |1953[ )’s early idea that “the death of a star in a supernova explosion may 
lead to the birth of a great number of new stars”, predicting that a ring of newborn 
stars would condense out of the expanding shell, it became a paradigm that one of 


the consequences of SNR/ISM interactions can be the formation of new stars. Van- 
|hala & Cameron|(|1998|l carried out numerical studies, concluding that shock waves 


with velocities in the range of 25 — 40 km s“ are capable of triggering collapse 
in molecular cores immersed in clouds interacting with SNRs. Melioli et al. ( |2006[ ) 
set constraints on the SNR radius and molecular gas density for which star forma¬ 
tion is allowed. If the SNR-cloud interaction is too strong, the cloud is completely 
destroyed; if it is too weak, the molecular core never collapses. In summary, the the¬ 
oretical conclusion is that if favourable conditions are given, triggered star formation 
can occur as a consequence of SNR/ISM interaction. 

Wootten ( 1978|l conducted one of the very first searches of star formation trig¬ 


gered by a SNR. From the analysis of an infrared object in a cloud near the SNR 
W44, the author concluded that it was an evidence of star formation estimulated by 
the expansion of the SNR. Since then, several searches have been carried out, in¬ 
cluding one based on multi-wavelegth data towards W44, where Paron et al. (20091 
demonstrated that the star formation was related to the nearby HII region G034.8- 
0.7 and not to the SNR. Other SNRs interacting with molecular clouds with indica¬ 
tors of active star formation have also been investigated, including W30 (|Ojeda-May 


et al.||2002|i, G54.1-H0.3 (|Koo et al.||2008|l, G357.7 -H0.3 ( [Phillips et al.|2009|l, IC443 


Odenwald & Shivanandan||1985[ et al.||2011| l, G18.8H-0.3 ( Dubner et al.pOO^ 
^aron et ar]|2012| l, etc. In all cases, even when the location of protostellar objects 
and young stellar object (YSO) candidates immersed in shocked molecular clouds 
were very promising, after comparing the characteristic timescale of star formation 
with the age of the SNRs, the general conclusion is that the stellar formation started 
before the SN explosion. 


Desai et al. (20101 carried out the most complete and homogeneous search for 


star formation related to SNRs by examining the presence of YSOs and molecular 
clouds in the environs of 45 SNRs in the Large Magellanic Cloud. After a very de¬ 
tailed analysis, based on different arguments (positions, densities, timescales, etc.) 
the conclusion is that there is none evidence of SNR-triggered star formation in the 
LMC. 


In conclusion, star formation is frequently seen near supernova remnants, but such 
physical association does not necessarily imply a causal relationship. Massive stars 
that end their lives exploding as SNe are formed in clusters or OB associations where 
the formation of new stars may continue and propagate outward for a prolonged pe¬ 
riod of time. It is then not surprising that core-collapse SNe are near young stars in 
star-forming environments. But apparently the SNR shocks result in an increment of 
turbulence that is not compatible with star formation. In fact, the agitation may be so 
violent that disperses the material, hindering further star-forming activity. This is an 
important open field of research that demands further theoretical and observational 
studies. 
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8 Comparison of radio emission with emission in other spectral ranges 


The emission in different spectral regimes traces material with different physical con¬ 
ditions. Briefly, optical filaments observed in most SNRs arise from shocked interstel¬ 
lar medium that is cooling radiatively, while in a few remnants the optical emision 
includes oxygen rich filaments, fragments of nearly pure ejecta launched from the 
core of the progenitor star during its explosion, or can be dominated by Balmer lines 
when a fast shock (velocities higher than ^ 200 km s“^) enters into partly neutral 
interstellar gas. The thermal X-ray emission contains essential information regarding 
the temperature, composition, distribution and ionization state of the material syn¬ 
thesized and ejected in SN explosions, and of the ambient matter swept-up by the 
supernova shock. If the origin of the X-ray emission is synchrotron, then its study 
can set powerful constraints on the role of the SN shocks in the production of cosmic 
rays. Infrared emission mainly marks the location of shock-heated dust. In summary, 
the body of multispectral observations is not only useful to understanding properties 
of the SNRs and their precursors, but also to explore the interstellar gas. 

The search for correlation between radio and IR emissions is a powerful tool for 
distinguishing thermal from non-thermal emission. The combination of the VLA sur¬ 
vey of the Galactic Plane at 330 MHz with the Midcourse Space Experiment (MSX) 
at 8/rm, served to discover 31 new SNRs and 4 candidates in the inner Galactic plane 
region, where the diffuse synchrotron emission and thermal HII regions cause more 
confusion. The SNRs are anti-correlated with dust emission, while HII regions are 
invariably surrounded by a shell of bright 8/rm emission ( Brogan et al.|2006] l. Also, 
this methodology has been applied to provide accurate estimate of the integrated radio 
flux density in remnants located in complex regions of our Galaxy where the contam¬ 
ination with thermal structures is high, as for example in the SNRs RX11713.7-3946 
( Acero et al.|2009| l and G338.3-0.0 ( Castelletti et al.|2011b| l. besides, the compari¬ 
son with radio emission from SNRs has been one of the prime methods to identify 
the nature of several 7— ray sources detected in the GeV and TeV ranges (more than 
twenty Galactic TeV sources are related with SNRs). 

Nowadays, the study of almost all SNRs is tackled from a multi-wavelength ap¬ 
proach. In what follow, we describe some sources in which the comparison of radio 
emission with emission in other spectral windows unveiled interesting results. 

G344.7-0.1 This SNR is a clear example where the comparison between radio, 
IR at 24^m and X-ray images elucidated the nature of a radio nebula located near 
the center of the remnant (Fig. [T^ . From a radio spectral study carried out in this 
region of the remnant using the Very Farge Array (VFA, NRAO) and the Australia 
Telescope Compact Array (ATCA) data at 1.4 GHz and 5 GHz, [Giacani et aL] ( |2011[ ) 
determined a mean radio spectral index of a ~ 0.3 for the nebula, a value compatible 
with those of radio PWNe. The combination of the emission in the three mentioned 
spectral ranges allowed the authors to rule out a PWN origin, concluding that the 
central bright radio feature is probably the result of strong shocks interacting with 
dense material. This interaction enhanced the infrared emission from shocked dust 
and favored particle acceleration, resulting in a flatter radio spectrum. 

G320.4-1.2 (MSH 15-52) A complex interacting system that has been observed 
throughout the electromagnetic spectrum is formed by the SNR G320.4-1.2, the en- 
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Fig. 12 Radio image of G344.7-0.1 at 1.4 GHz taken from 
the location of the radio nebula near the geometric center of the remnant 


'2011 


. The green square shows 


ergetic pulsar PSR B1509-58 and its PWN, and the Hll region RCW 89 (Fig.[T^. 
SNR G320.4-1.2 has an unusual radio appearance consisting of two distinct compo¬ 
nents: towards the northwest, a bright centrally concentrated source with a ring of 
radio clumps in coincidence with the optical nebula RCW 89, and towards the south¬ 
east a fainter partial shell (Fig. [T3]Le/fj. The pulsar PSR B1509-58, detected in radio. 
X-rays and in 7— ray is located near the center of G320.4-01.2. In the X-ray domain, 
this system has a complex picture with a different morphology (Fig. [l^Rig/if). The 
emission is dominated by a jet-like structure emerging from the pulsar, the PWN, 
and extended emission, thermal in origin coincident with the northern radio compo¬ 
nent of the SNR and with the Fill region RCW 89. Gaensler et al. ( 1999| l proposed 
that this complex system can be explained as the result of a low-mass or high-energy 
explosion occurring near one edge of an elongated low density cavity, which was 
confirmed by observations of the neutral gas in the region by Dubner et al. (2002 1 . 
The pulsar appears to be generating twin collimated outflows, the northern part of 
which interacts with the SNR producing the collection of the clumps observed in ra¬ 
dio and in the X-ray band. Gaensler et al. (20021 later reported the discovery of the 
long-sought PWN on the basis of comparison of new Chandra X-ray data with radio 
morphology. 

G292.0+1.8 This remnant is another interesting source that has been subject to 
multi-wavelength observational campaigns. It is a textbook example because it has 
all the properties of a core-collapse SNR: a central pulsar (PSR J1124-5916) feed- 
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Fig. 13 Left: Radio continuum image of the SNR G320.4-1.2 at 1.4 GHz taken from 
Right: X-ray image by Chandra jGaensler et al.|200^. Credit Image NASA/MIT/B. Gaensler et al. 


2002 


ing a PWN, metal-rich ejecta, shocked CSM, and a blast wave. For a summary of 
previous observations, see Ghavamian et al. ( |2012| l. Briefly, it is characterized by a 
remarkable network of filaments seen in optical and in X-rays, which are fragments 
of supernova ejecta ( [Winkler & Long|2006[[Park et al.|2002[|2004]|2007| l. The X-ray 
image also reveals bright filaments running from east to west, named the “equatorial 
belt”, probably shocks propagating in circumstellar material ( Park et al.|2002 1 . In the 
radio band, a mutifrequency study of the remnant performed with ATCA at 20, 13 
and 6 cm by Gaensler & Wallace ( 2003| l, revealed the presence of a polarized and 
flat-spectrum bright central core, representing the radio PWN, surrounded by a circu¬ 
lar fainter plateau with steep-spectrum which represents the SNR shell. Curiously, a 
series of radial filaments with much flatter spectrum are located over the plateau. The 
comparison with the X-ray emission allowed Gaensler & Wallace ( |2003| ) to propose 
that Rayleigh-Taylor instabilities near the SNR contact discontinuity originated the 
formation of such filaments (Fig.3 in Gaensler & Wallace|2003 1. 


Gl.9+0.3 It is the youngest remnant detected in the Milky Way (age ~150 yr), 
opening as such the opportunity to study a Galactic SNR on its very early develop¬ 
ment. The morphology of G 1.94-0.3 can be described as an almost complete shell in 
radio and in X-rays ( Green et al.||2008] Reynolds et al.||20081 l, but the comparison 
between the distribution of the emission in both spectral ranges revealed a peculiar 
anti-correlation. While the radio remnant is clearly brighter along the northern bor¬ 
der, the X-rays has a notable bilateral east-west symmetry including two extensions 
(“ears”) not detected in the radio band (Fig.[l4jl. It is one of the few shell SNRs with 
an X-ray spectrum dominated by synchrotron emission ( Reynolds et al.|2009] l. But 
this is not the only origin of the X-ray emission, Borkowski et al. ( 2010| l reported the 
presence of lines of Fe, Si, and S in small regions of the northern limb with spec¬ 
troscopic velocities of about 14000 km s“^ and a line at 4.1 keV, identified as due 
to ^^Sc (the first firm detection of this element in a SNR). Based on the presence 
of Fe lines, high velocities, absence of PWN, and bilaterally symmetric non-thermal 
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X-ray emission as in SN 1006,|Borkowski et al.|P0131l suggest a Type la origin for 
G1.9H-0.3. 



Fig. 14 Chandra image of Gl.9+0.3 (in green) jBorkowski et al. 


2013 


superimposed with VLA radio 


image at 4.9 GHz (in red) as reprocessed from archival VLA data by the authors of this paper 


9 Supernova Remnants in the Magellanic Clouds 


The Magellanic Clouds are an ideal laboratory to investigate the properties, both 
global and individual, of SNRs across the whole electromagnetic spectrum. The fact 
that the distances are known (49 kpc for the LMC and 59 kpc for the SMC), that 
the LMC is nearly face-on, so that all remnants are nearly co-distant, and that even 
for the rounder SMC there is a range of only ± 5 kpc with respect to the mean 
distance, allows for accurate estimation of physical parameters of the remnants. At 
the same time, the proximity permits radio observations with good angular resolution 
revealing the structure of individual objects and facilitating accurate multi-spectral 
comparisons. 

One of the very first identifications of a SNR in the Magellanic Clouds was re¬ 
ported by Mathewson & Healey ( 1963) 1 based on radio observations carried out with 
the 210-ft. steerable reflector at the Australian National Radio Astronomy Observa¬ 
tory. It was followed by several searches using the Molonglo Observatory Synthesis 
Telescope (MOST) that permitted the identification of new members based on the 


combination of radio observations with optical and X-ray data. Mathewson et al. 
( 1983| l compiled a catalogue of SNRs in the Magellanic Clouds with 25 confirmed 
remnants in the LMC and 6 in the SMC. Later on, the search and studies of SNRs 
in our neighbouring dwarf galaxies took advantage of the Australia Compact Tele¬ 
scope Array (ATCA) that improved the angular resolution by factors between 5 and 
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35 compared with the best previously available radio data. Several works reported 
discoveries of new members and detailed studies of already known SNRs (e.g.jAmy 


|& Ball|[T993l [Manchester eTan|1993a|b[ [Dickel et al.|[T993l jPickel & Milne||1995 

[T998letc.) 


Badenes et al. (2010|l compiled from the existing literature a catalogue of multi¬ 


wavelength observations with 77 confirmed SNRs in the Magellanic Clouds (54 in the 
LMC and 23 in the SMC), arguing that this list comprises a fairly complete record of 
SNe that exploded over the last ~ 20 kyr. After Badenes published the catalogue, at 
least seven new cases were confirmed in the LMC (jGrondin et al.|201^ Maggi et al. 


2012[ [Kavanagh et al.|2013| [Bozzetto et al.|2013| [Maggi et al.|2014 1. Of all the 84 

SNRs identified to date, only 4 of them (SNR J005.9-7310 in the SMC, and MC SNR 
J058-6830, MC SNR J0511-6759 and MC SNR J0517-6759 in the LMC) have not 
been detected in the radio band, possibly because they are old and faint ([Maggi et al.| 


[M4l ). 

Badenes et al.[ ( [2010| l carried out statistical studies of the SNRs in the Magellanic 
Clouds concluding that the size distribution of remnants is approximately flat with a 
cutoff at r^ 30 pc. The authors propose that most of the SNRs are in the Sedov phase 
of evolution, quickly fading below detection as soon as they reach the radiative stage. 

Among all SNRs in the Magellanic CLoud, the radio remnant of SN 1987A de¬ 
serves special attention because it offered the possibility of tracing in great detail its 
temporal evolution. In what follows, the main aspects of this source are summarized. 


The radio remnant of the SN 1987A SN 1987A was the first naked-eye SNe event 
since the invention of the telescope. It exploded in 1987 February 23 in the Large 
Magellanic Cloud. Briefly, the progenitor of SN 1987A, Sk-69°202 with an initial 
mass estimated in ~ 20 M©, is believed to have evolved from a red supergia nt (RSG) 
into a blue supergiant (BSG) approximately 2 x 10^ yr prior to the explosion (Crotts & 


Heathcote|2000|l. Optical imaging revealed a complex circumstellar medium (CSM) 


surrounding the explosion, consisting of a triple-ringed structure (Burrows et al. 


1995|and reference therein). The two outer rings appear to be the cap of an hourglass 


shaped structure enveloping the SN itself. The inner ring, also referred to as the equa¬ 
torial ring, is believed to represent an equatorial density enhancement in the CSM, 
located at the interface between a dense wind emitted from an earlier red-giant phase 
of the progenitor star and a faster wind emitted by the star in more recent times. 

In the radio band, the emission from the SN 1987A was detected at 843 MHz 
with the MOST two days after the SN event. This emission reached its maximum 
value of around 140 mJy 4 days after the explosion and then decayed rapidly to 
become undetectable less than a year later ( [Turtle et al.[|1987) . This radio outburst 
has been explained as the consequence of the rapidly moving, low density BSG wind, 
which produced only a short-lived period of radio emission when hit by the SN shock 
( Storey & Manchester|1987]|Chevalier|1998| l. 

Approximately 1200 days after the explosion, radio synchrotron emission was 
again detected, in this case by both ATCA and MOST ([Staveley-Smith et al.|1992 


Ball et al.|2001 respectively), marking the birth of the radio remnant. From the first 


radio detection, the radio emission has been regularly monitored every 1-2 months at 
1.4, 2.4, 4.8 and 8.6 GHz using various array configurations of ATCA ([Manchester[ 
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|et al.||2002{ |Staveley-Smith et al.||2007[ [Zanardo et al.||20l0] and reference therein). 
At all these frequencies, the flux density has been steadily increasing over time and 
from ^ day 3000 after the explosion is growing exponentially, which is attributed 
to increasing efficient particle acceleration processes ( |Zanardo et al.|[2OT0 1 . On the 
other hand, the highest angular resolution observations (~ O'M) of the SNR obtained 
with the Australian Large Baseline Array (LBA) at 1.4 and 1.7 GHz revealed two 
extended lobes with an overall morphology in good agreement with that at lower 
angular resolutions. In addition, small-scale structures were found in the brightest 
regions in both lobes (|Tingay et al.|2009t Ng et al.|201 11. 


Recently, Ng et al. ( 2013| l reported on the study of the evolution of the radio 
morphology of the remnant of the SN 1987A covering the period January 1992 - May 
2013 (day 9568 after the explosion). The data were acquired with ATCA at 9 GHz 
with an angular resolution of 0".4. The remnant presents a double lobe ring, with 
an asymmetric surface brightness distribution being the eastern lobe brighter than 
its western counterpart (Figure p3] l. From the analysis of this database the authors 
point out that from day 7000 the asymmetry began to decline, such that the overall 
geometry is evolving towards a ring structure, suggesting that the remnant has entered 
in a new evolutionary stage where the forward shock has fully engulfed the entire 
inner ring and is now interacting with the densest part of the circumstellar medium. 


Observations at higher frequencies (18, 36, 44 and 94 GHz) undertaken with 
ATCA ( [Manchester et al.| | 2005 [ [Potter et al.| | 2009 | [Lakicevic et al.[ | 2012[ [ Zanard o[ 
et al.pOll I, followed by observations at 110, 215, and 345 GHz with the Atacama 


Millimeter/Submillimeter Array (ALMA), show the same asymmetric double lobe 
morphology (Indebetouw et al. 2014 Zanardo et al. 2014| l. Also ALMA observa¬ 
tions ( jlndebetouw et al.|2014 1 revealed the presence of one of the largest masses of 
cold dust (greater than 0.2 M©) ever measured in a SNR. The dust emission is con¬ 
centrated at the center of the remnant, and has probably formed in the inner ejecta. 


The evolution of the radio spectral index has been tracked over the years. At 
present the mean spectral index estimated between the frequencies 18 and 44 GHz is 
a = 0.8, with values in the range 1.1 to 0.3 across the source. The brightest region on 
the eastern lobe is slightly steeper than the mean value, while flatter values, between 
0.5 and 0.3, are found toward the central and the center-north regions of the SNR 
( Zanardo et al.|20r3| l. The spectral index distribution of the synchrotron component 
was mapped in the frequencies range 44 to 672 GHz with data obtained with ATCA 
and ALMA instruments ( [Zanardo et al.|2014l ). From radio to far infrared, the authors 
estimate a mean spectral index for the synchrotron main component of 0.73, while 
a value ~ 0 was estimated around the central region and between 0.1 and 0.4 in the 
western half of the remnant. Zanardo et al. ( 2014| l conclude that central region might 
be a PWN in the interior of the remnant, while the steepening of the spectral indices 
over the eastern lobe might be indication of a local spectral break. 
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10-1 
RA offset (arcaec) 


Fig. 15 9 GHz images of SN 1987A in the interval 1992-2013 obtained with a FWHM of 0".4 (figure 
provided by|Ng et al.|2013) 


10 Conclusions and prospects of studies with the new generation 
radiotelescopes 

Almost seven decades after the first detection of the radio emission coming from 
SNRs, great advances in the knowledge of their properties have been achieved. De¬ 
tailed studies of the radio emission associated with SNRs have proven to be an excel¬ 
lent tool to recognize morphological characteristics, delineate the location of shock 
fronts, identify sites of interaction with dense clouds of the ISM, locate places of 
particle acceleration, and discover PWNe even in cases where the radio pulsar is not 
detected and the nebula is the only trace of its existence. They are also powerful 
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probes of the magnetic field, providing information on the orientation and sites of 
higher compression. Furthermore, the possibility of resolving individual radio struc¬ 
tures down to spatial scales of arcsec or even sub-arcsec, like in the expansion studies 
performed in Tycho’s and SN1006 SNRs or the current studies of the Crab Nebula, 
is a powerful resource to investigate the dynamical evolution of SNRs. 

Over the past years, key observations performed with single-dish telescopes and 
radio interferometers like ATCA, DRAG, VLA, VLBA, GMRT, etc. provided re¬ 
markable insights on radio SNRs. However, there still persist controversial issues and 
unsolved questions requiring further research from observational as well as from the¬ 
oretical approaches. In brief, particular issues deserving further observational studies 
are: the deficiencies of the spectral database, the generally poor quality of spatially 
resolved spectral indices over wide frequency ranges, the scarcity of well determined 
distances, the incomplete knowledge of extended old SNRs, the missing (young and 
old) Galactic remnants, the consequences at all levels of the shock/cloud interactions, 
the alignement and strength of the magnetic field, and where and how particles are 
accelerated. From the theoretical point of view, the comprehension of SNRs would 
benefit from progress in the research of mechanisms of particle acceleration capa¬ 
ble of explaining the variety of observed spectral indices, as well as from detailed 
modelling of the dynamical evolution of SNRs expanding in inhomogeneous envi¬ 
ronments. 

Other motivation to increase the quality and accuracy of radio measurements is 
to match the ongoing revolution in high energy (X-ray and 7— ray) studies. It is well 
known that recent surveys of the Milky Way with space and ground-based 7— ray de¬ 
tectors revealed hundreds of high-energy and tens of very-high-energy 7— ray emit¬ 
ters representing several Galactic source populations, including numerous SNRs and 
PWNe whose emission mechanisms remain poorly constrained. The fundamental na¬ 
ture of many of the high energy sources could greatly benefit from radio observationes 
that are only recently becoming technically feasible, e. g. with the JVLA or ALMA. 

Another area with pending studies is the investigation of faint SNRs, where the 
lack of precise information from radio data may affect the comprehension of high- 
energy phenomena and, at large, the question of the contribution of SNRs to the 
overall flux of Galactic cosmic rays. This is, for example, the case of the SNRs RX 
J1713.7-3946 and RX J0852.0-4622, two very interesting shell-type X-ray SNRs with 
bright 7— ray emission (at 1 TeV RX J0852.0-4622 is as bright as the Crab Neb¬ 
ula), but with poorly constrained radio parameters. These sources are radio weak, 
confused on parts with thermal emission and with significant Galactic background 
variation, characteristics that all together are detrimental to the accurate estimate of 
radio parameters. For cases like these, the new low-frequency synthesis arrays like, 
low-frequency array (LOFAR) in Europe, the Murchison widefield array (MWA) in 
Australia, and the Low-Band JVLA system in the USA, are important new tools. 
Their good sensitivity and angular resolution is key to imaging and differentiating the 
non-thermal radiation (bright at low radio frequencies) from optically thick and thin 
thermal emission from HII regions (which are intense at high-frequencies). Low radio 
frequency studies are also indispensable to improving the precision of the continuum 
spectra of all SNRs. Accurate low-frequency data anchor the spectrum, avoiding the 
uncertainties of extrapolation from higher frequency data which render them insensi- 
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live to curvature, an important piece of information to test the predictions of particle 
acceleration theories. In addition, sensitive high-resolution meter-wavelength studies 
of selected Galactic regions, have already proven to be extremely useful to discover 
new SNRs. The detection of many new SNRs would alleviate the incompleteness of 
the current Galactic SNRs census with important implications for the SNe rate and 
energy input into the ISM. A deep meter-wavelength census of Galactic SNRs with 
the JVLA is currently under way, as are LOFAR studies of selected SNR complexes. 

In the other extreme of radio frequencies, observations of radio-continuum emis¬ 
sion of bright SNRs in the tens and hundreds of GHz regime using ATCA, ALMA, 
and JVLA, among other facilities, have come to fill in the gap between radio and IR, 
offering significantly improved constraints on SNRs spectral energy distributions. 

The new aperture-synthesis telescopes, either in operation, in commissioning, or 
in construction phase (EVLA, SKA and its precursors ASKAP in Australia, MeerKAT 
in South Africa, LWA in the USA, and FAST in China, plus the already mentioned 
MWA) will also bring important advances in understanding the properties of the mag¬ 
netic field in SNRs. They will provide high angular resolution and will be able to 
produce wideband multichannel data amenable to rotation measure (RM) synthesis. 
Wide area polarization surveys will be very helpful to establish the presence of mag¬ 
netic field, its direction, strength and spatial scale in a large sample of Galactic SNRs 
and in the Magellanic Clouds. In addition, very long baseline interferometers, like 
the European Network eEVN, are the appropriate instruments to resolve the structure 
and investigate the expansion of young SNRs located in nearby galaxies other than 


the Magellanic Clouds, following the pioneering studies conducted by Pedlar et al. 


( |1999| l in M 82 and |Bartel et al.| ( [T994| ) and |Marcaide et al.] ( |1995| ) in M 81, and the 
more recent from |Bietenholz et al. ( 2010a| l in NGC 4449, Bietenholz et al.| ( f2010b[ ) 
in NGC 891, etc. 

Another chapter connected with the investigation of radio SNRs that deserves 
special attention is the study of the atomic and molecular emission from the surround¬ 
ing ISM. Such studies have become more and more important to understanding the 
properties of the SNRs across the whole electromagnetic spectrum. Large-scale HI 
surveys, like SGPS (Southern Galactic Plane Survey) in the southern hemisphere and 
VGPS (VLA Galactic Plane Survey) in the north, have been very helpful to explore 
in emission and absorption the characteristics of the ISM towards a large number of 
SNRs close to the Galactic plane. Also intermediate resolution CO surveys, like the 
GRS (Galactic Ring Survey of ^^CO in the northern sky) and the ^^CO NANTEN 
4-m dish observations in the southern sky, have been important to localize possible 
sites of SNR/MC interaction. For sources at higher Galactic latitudes, the atomic and 
molecular studies require dedicated observations. The necessity, however, of more 
detailed studies of the surrounding ISM for a large number of SNRs is evident. An 
ample, good-quality atomic database will allow us to refine distance estimates and, 
particularly, HI data acquired with very high angular resolution and sensitivity can 
be the way to detect the HI shell predicted to form behind the expanding shock front 
when radiative cooling is sufficient to recombine a detectable amount of atomic ma¬ 
terial. On the other hand, molecular data at sub-arcmin angular resolution, as can be 
obtained with single-dish of intermediate size, like the 12-m antenna of the future 
Argentina-Brazil radiotelescope Long Latin American Millimetre Array (LLAMA), 
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will be an adequate step forward following the NANTEN southern sky survey. The 
data will have the sensitivity and angular resolution suitable to cover large areas iden¬ 
tifying compact molecular clumps that might have been overrun by SNRs. 

In view of the present radio astronomy landscape, with new instruments already in 
use or coming into operation in the next future, plus upgrading of existing instruments 
with the latest technological developments, the panorama of future research of radio 
SNRs is very promising. 
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